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ABSTRACT 


In Active Galactic Nuclei (AGN) the luminosity is so intense that the effect of 
radiation pressure on a particle may exceed the gravitational attraction. It was shown 
that when such luminosities are reached, relatively cold (not completely ionized) 
thermal matter clouds may form in the central engines of AGN, where most of the 
luminosity originates. 

We show that the spectrum of emission from cold clouds embedded in hot rela- 
tivistic matter is similar to the observed spectrum. We also show that within the hot 
relativistic matter, cold matter moves faster than the speed of sound or the Alfven 
speed, and shocks form. The shocks provide a mechanism by which a localized per- 
turbation can propagate throughout the central engine. The shocked matter can emit 
the observed luminosity, and can explain the flux and spectral variability. It may also 
provide an efficient mechanism for the outward transfer of angular momentum and 
the outward flow of winds. 

With observations from X-ray satellites, emission features from the cold and hot 
matter may be revealed. Our analysis of X-ray data from the Seyfert 1 galaxy MCG - 
6-30-15 over five years using detectors on the Ginga and Rosat satellites, revealed some 
interesting variable features. A source with hot matter emits non-thermal radiation 
which is Compton reflected from cold matter and then absorbed by warm (partially 
ionized) absorbing matter in the first model, which can be fit to the data if both the 
cold and warm absorbers are near the central engine. An alternative model in which 
the emission from the hot matter is partially covered by very warm matter (in which 
all elements except Iron are mostly ionized) is also successful. In this model the cold 
and warm matter may be at distances of up to 100 times the size of the central engine, 
well within the region where broad optical lines are produced. The flux variability 
is more naturally explained by the second model. Our results support the existence 
of cold matter in, or near, the central engine of MCG -6-30-15. Cold matter in the 
central engine, and evidence of the effects of shocks, is probably forthcoming with 
future X-ray satellites. 
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Chapter 1 

INTRODUCTION 

A 

In this chapter we introduce the current theories and observations of the central 
engines of AGN. In the first section partially successful models of the central engines 
of AGN are presented. The second section is devoted to observations of AGN in the 
X-ray waveband, in which most of the radiation from the central engines is probably 
emitted. In the third section some of the problems in theoretically modeling the 
central engines of AGN are presented, possible solutions suggested in this thesis are 
outlined and our analysis of X-ray data from an AGN is discussed. 

1.1 INTRODUCTION TO THEORY OF AGN 

1.1.1 What are AGN? 

} 

Active galactic nuclei (AGN), which are the extraordinarily luminous centers of 
some galaxies, axe among the most enigmatic phenomena in astrophysics. AGN may 


j 



be as small as the solar system, but have energy outputs which, in some cases, rival 
the luminosity of the entire host galaxy. Approximately 10% of all galaxies are AGN, 
and at least 50% of galaxies show some evidence of activity in their nucleus. As of 
this date there is no one cohesive explanation to the source of this phenomenon that 
can be directly corroborated by observations. AGN are therefore characterized by 
observations. 

In table 1 some of the observable properties of AGN are presented. Only properties 
related to the galactic nucleus are included. Other basic properties can be found 
in Woltjer 1990, Netzer 1990, Mushotzky, Done and Pounds 1993 and references 
therein. A partial list of AGN and host galaxies types is presented in the first two 
columns: Seyferts are optically resolved spiral galaxies with an AGN. Seyferts of 
type 1 have both broad and narrow emission lines in the infra-red (IR) - ultraviolet 
(UV) wavebands, and Seyfert 2s have only narrow emission lines. Seyferts comprise 
about 5% of all galaxies. AGN which could not be optically resolved when they 
were first observed are usually referred to as quasars. (In the past few years some 
quasars were optically resolved and shown to be AGN.) The radio quiet quasars are 
called quasi-stellar objects (QSO). In optically resolved elliptical galaxies AGN are 
usually radio galaxies and BL LAC objects (whose spectrum resembles the spectrum 



of strong radio galaxies, but lack emission lines and are highly variable). The radio 
loud quasars are quasi stellar radio sources (QSR) and quasars which exhibit strong 
optical variability are called optically violent variables (OVV, mostly a subgroup of 
QSR). Increased activity in the nuclear region may also be in the form of increased 
rate of star formation (in star burst galaxies), strong emission lines (in LINERS and 
nuclear H II regions) and increased IR emission (strong IRAS galaxies). 

AGN morphology is based primarily on luminosity and spectrum, although the 
property which distinguishes AGN from most other astrophysical objects is their high 
luminosity per volume. We therefore included in the next four columns of table 1 
some information on observed properties which are directly related to the luminosity 
and size. The estimated total luminosity in column 3, L tot , includes emission in all 
wavelengths, whereas in column 4 we only include Lx, the emission in the dominant 
x-ray waveband. In column 5 the estimated doubling time scale A t (the minimum 
time the X-ray flux takes to double in magnitude) is presented. Using columns 4 
and 5 the compactness of the source can be defined as / = Lx<tt/{ Atm e c 4 ) (Guilbert 
Fabian and Rees 1983), where <tt is the Thompson cross section, m e is the electron 
mass and R = cAt is usually referred to as the size of the central engine (see sections 
§1.3 for a more precise physical definition of the central engine). The compactness 
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parameter is only appropriate for the first three categories of AGN. In Seyfert 2s the 
nucleus is most probably obscured, and the compactness parameter is meaningless. 
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Table 1. AGN: Observations That Can be Related to Central Engines 


AGN type 

Host Galaxy 4 

r b 

l^tot 

Lx c 

At d 

l e 


Seyfert Is 

spiral 

lQ‘12-45 

10 4 1 —44 

1Q3-5 

io~-~ 

10 3 

Sevfert 2s 

spiral 

10 42- ' 15 

10-10—12 

~ 10 7f . 

o 


QSO 

spiral(’) 

1Q44-4T 

1Q43-4C 

10 5-T 

1CT 2 - 

I0 3 

Radio 

elliptical 

~ 10‘--“ 5 s 

S 

10-1-7? 

? 


Strong Radio 

elliptical 

~ 10‘ l2_ ‘ u * 

S 

•? 

? 


QSR 

elliptical(?) 

~ 10 43- ‘ ,T * 

& 

? 

? 


BL Lac 

elliptical 

h 

h 

~ 10“ 

? 


ovv 

all? 

h 

h 

10 4 

? 


LINER 

spiral 

i 

i 

? 

’> 


HII regions 

spiral' 

i 

? 

? 



IRAS 


~ 10‘ 15j 

? 

? 

? 



a QSO and QSR are mostly unresolved. 

b in erg sec -1 , estimated for the galactic nuclei. 

c in erg sec -1 in the 2 — lOKev range, estimated for the galactic nuclei. 
d in seconds. Upper limit is an estimate. 

'Compactness of central engine, lower limit: l ~ L z <Tt/ R m e c 3 

f Some evidence for variability between x-ray observations 

®Most of the observed flux is not from the central engine 

k Most of the observed flux is probably from an observer oriented jet. 

'Galactic nuclei emission can only be estimated in LINER and HII regions. 

j Most emission in IR. 
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1.1.2 The Importance of AGN 

Better understanding of AGN may help answer a number of questions. Among 
them: 

1. Galactic morphology: AGN may help us understand how galaxies form and evolve. 
A large fraction of all galaxies have at least some AGN features, which implies that 
most galaxies may have had AGN in their centers at some point in their history 
(Blandford 1990, and references therein). Emission from highly red shifted quasars 
imply the existence of nuclear galactic activity at an epoch in which galaxies were 
not even supposed to have formed, according to many scenarios of galaxy formation. 
Understanding AGN will therefore help in understanding galactic evolution. 

2. Cosmology: If we can understand AGN, quasars may become the perfect ‘stan- 
dard candles’, which can be used to determine the geometry of the universe . Us- 
ing a standard candle of known luminosity Lq at a ‘luminosity distance’ = 
(Lq/AtcF) 1 / 2 , the cosmological parameters can be determined by observations us- 
ing H 0 di = z + 0.5(1 — qo)z 2 + ..., the approximate solution of Einstein equations 
for an expanding matter dominated universe. Here F and z are the observed flux 
and redshift, respectively. The Hubble constant Ho and deceleration parameter qo 
may be found by using quasars, the highest redshift objects in the universe. These 



parameters, if determined, may reveal the geometry of the universe (Woltjer 1990, 
Schutz 1985). 

3. High energy physics: AGN may serve as high energy labs with which one can 
enhance our understanding of high energy phenomena. There are reasons to believe 
that very efficient particle acceleration mechanisms are present in AGN centers, and 
that plasmas with temperatures unaccessible in labs are present near their centers 
(Blandford 1990, Blandford and Eichler 1987). 

4. General relativity: Understanding AGN may allow us a look into the very nature 
of the general theory of relativity. Some effects that may be observed in the future 
are the Lense-Thirring precession, energy extracted from a rotating black hole or, 
possibly, gravitational radiation from binary black holes in the central engine of an 
AGN (Blandford 1990 and references therein). These effects, when observed, will 
consolidate our understanding of general relativity. 

5. The large scale structure of the universe: The spectrum and isotropy of the X- 
ray background are probably the signature of a period less than a billion years after 
the big bang. In that period more galaxies had active nuclei. It was shown that 
AGN may account for most of the X-ray background. An early evidence of the large 
scale structure may therefore be available (Fabian and Barcons 1992, and references 



therein). 

In order to advance towards a solution to these questions better understanding 
of AGN must be developed. In particular, emission of light and particles from the 
vicinity of the AGN center must be understood. These processes can only be explained 
in the context of a model for the central driving force of an AGN, the central engine. 

1.1.3 Possible Candidates for the Central Engines of AGN 

It is now generally accepted that super massive compact objects lurk at the centers 
of AGN (see Blandford 1990, appendix 1). Gravitational potential energy of matter 
which falls onto these objects is released in the form of radiation. Although the 
gravitational behavior of the massive objects is yet unknown, and the mechanisms 
by which these objects generate the incredible observed phenomena are not yet fully 
explained, the matter accreting super massive compact object model is the most 
successful in describing the phenomena. 

Traditional candidates for the AGN central engines, such as dense clusters of 
normal luminous stars and clusters of neutron stars are probably rejected by the 
recent observations with the Hubble Space Telescope (HST) of a single central super 
massive object (see Ford et al. 1994). The Chandrasekar limit on masses of stars 
makes it impossible to construct a super massive compact object made of baryonic 
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matter. Another alternative model is the star-burst model (Terlevich et al. 1992, and 
references therein), in which the energy is generated by violent star formation activity 
in the innermost regions of AGN. Although the spectra of AGN may be explained in 
this manner the variability is almost impossible to explain (Green 1993, and references 
therein). Its also hard to reconcile star-burst models with HST observations. 

The compact objects are assumed to be super massive black holes (SMBH) in 
the so called ‘best bet’ model for AGN (Salpeter 1964, Lynden-Bell 1969, Rees 1984, 
Blandford 1990). We hereafter refer to the central compact objects as SMBH, as 
is often done in the professional literature. SMBH are the central objects in the 
models described in the following sections. (But note the quote from Blandford 
1990: “...various observational discoveries are increasingly hard to interpret in terms 
of alternative models... however... the evidence [for black holes] is weaker than for 
stellar mass black holes in X-ray binaries...) 

Black holes are vacuum solutions of Einstein equations with an event horizon, from 
which light and particle can not escape (Wald 1984). Schwarzschild and Kerr family of 
solutions span all possible solutions for vacuum stationary axisymmetric spacetimes 
with a regular event horizon (Misner, Thorne and Wheeler 1972, Robinson 1975, 
Mazur 1982). The spherically symmetric Schwarzschild solution is a special case of 
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the axisymmetric Kerr solution. A positive change in the mass of the central compact 
object, AM > 0, is required in both of the Schwarzschild and Kerr solutions, if the 
efficiency of the accretion process is less than 100 %. 

Black holes are more efficient in conversion of the mass to energy than other as- 
trophysical processes. For example, nuclear fusion of hydrogen to helium can convert 
up to 0.7% of the rest mass to energy. Near Kerr SMBH up to 40% of the rest mass 
of the gas accreted onto the black hole can be converted into energy. 

The release of gravitational energy from matter that falls onto the SMBH is via 
viscous processes in either accretion disks or semi-spherically accreted matter (see 
Frank, King and Raine 1992; Guilbert and Rees 1987; chapters §2 and §5 of this 
thesis). In accretion onto black holes most of the gravitational energy is released 
within the inner 10/? p where R g = GM/c 2 and M is the mass of the SMBH. This 
can easily be seen if the Newtonian potential energy becomes kinetic energy which 
is radiated away, such that GMmu /(6R g ) — GMmn/(lOR g ) ~ 0.06mj/c 2 , or 6% of 
the rest mass is converted into energy, where m// is the mass of hydrogen. The inner 
boundary is chosen to be the last stable orbit in a Schwarzschild geometry. Accurate 
calculations of the total energy which can be radiated by an accreted particle yield a 
similar number. The efficiency of the accretion process is not only higher than that 
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of hydrogen fusion, for example, but is also more stable than fusion when limited to a 
size 10/? 3 . We hereafter refer to the inner 10R 9 as the ‘the central engine’ as is often 
done in the professional literature. Accretion processes are assumed to be the main 
source of the energy emitted from the central engine. 

The limited supply of mass to be accreted onto the SMBH in the center of AGN 
poses some problems for the best bet model. For a thorough discussion of this problem 
the reader is referred to Rees 1990. 

Direct emission from the central engine may give us evidence of a SMBH. However, 
the spatial resolution of present day telescopes can not tell us whether most of the 
radiation originates in the central engine or elsewhere, even in the nearest AGN. 
Redshifted emission lines may indicate emission from the vicinity of a SMBH, if 
observed. Variability may be observed which could be explained by gravitational 
lensing modification of emission from matter that follows geodesics in the central 
engine (Rauch and Blandford 1994). 

A combination of observations of the spectrum, flux variability and spectral vari- 
ability can be used for reconstruction of the possible geometrical configuration of the 
central engine. The results of such observations are described in the next three sec- 
tions. A way to distinguish the central engine radiation from competing emissions 
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regions may be found in the near future. 

1.1.4 Spectra of Emission from the Central Engines of AGN 

The overall continuum spectra of most AGN can be characterized as power-law 
with bumps super-imposed in the IR and UV, or extreme UV (EUV) and X-ray 
wavebands. (The IR bump probably does not originate in the central engines of 
AGN. See Blandford 1990.) Emission and absorption features are also present at 
photon energies of up to ~ lOkeV. 

In figure 1-1 we show the spectrum in wavebands from radio to 7 rays of the quasar 
3C273. The abscissa has units of v and the ordinate has units of log i/F u , where v is 
the photon frequency and F„ is the flux per frequency. The reason for this choice of 
coordinates is that for a power-law spectrum, of the form F„ oc v~ T , the integrated 
flux between two frequencies, u\ and 1 / 2 , is F = CXog^jvi) when T = 1 , where C 
is a constant. The total emission in each photon frequency decade is therefore the 
same, and a ‘flat’ curve appears on a graph of the above format. It is evident that 
the AGN flux in figure 1-1 has an almost flat spectrum. Black body stellar emission 
of a galaxy appears as a reference for comparison. 
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Figure 1.1: The spectrum quasar 3C273. AGN emit similar integrated power in the IR 
7 -ray wavebands (see text). (Adapted from Green 1993.) 
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If a SMBH is at the center of the AGN then most of the radiation originates in the 
central engine, as demonstrated in section 3 and chapter 2. The radiative processes 
responsible for the radiation that emerges from the central engine probably include 
synchrotron radiation, Compton scattering, radiative reprocessing in pair plasmas, 
and may includes atomic radiative transitions and black body emission. 

Synchrotron radiation emitted from relativistic electrons in magnetic fields should 
affect the spectrum if the synchrotron cooling time is shorter than the typical accretion 
time. Synchrotron emission, which should result in radio polarization, was proven to 
be important in most AGN by measurements of polarization. Measurements of mag- 
netic fields in galactic centers yield values which require that synchrotron radiation 
be an important radiation mechanism on all scales. These magnetic fields which 
probably originate in stars which are being disrupted by the SMBH are increased as 
the accreting plasma becomes denser. Near the central engine the magnetic fields 
are strong enough to result in cooling time scales which are much shorter than the 
accretion time scales. 

Synchrotron radiation in astrophysics usually result in an almost flat power law 
spectrum, similar to the underlying spectrum of most AGN. The synchrotron radia- 
tion in astrophysics usually comes from a power-law distribution of electrons veloci- 
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ties, which may result from acceleration mechanisms such as the Fermi-shock process 
in shocks. The resultant spectrum is an almost flat power-law spectrum, which almost 
with a low frequency cutoff due to synchrotron self absorption. This self absorption 
may be associated with the radio-quietness of some compact AGN. 

The effects of Compton scattering on the spectrum are noticeable when the cooling 
time of photons of energies comparable to m e c 2 is shorter than the time in which a 
photon is crossing the emission region, t x = R/c. The spectrum may also be affected 
by Compton scattering of low energy photons by relativistic electrons when the time 
scale for that process is shorter than t x . Compton heating and cooling are important 
if the overall Thompson (or Klein Nishina) optical depth is large, which happens 
when the source is sufficiently compact (/ > 1, see table 1). Compton cooling, usually 
referred to as ‘reflection’, is important if the Thompson optical depth of cold matter 
in the photon’s path is larger than 1. The result of Compton cooling on a power-law 
spectrum is a decrease of the 7-ray photon flux and an increase in hard X-ray photon 
flux. Compton processes may, therefore, be responsible for the high energy cutoff and 
hard X-ray bump in the spectra of radio quiet AGN (Mushotzky, Done and Pounds 
1993 ). 

Pair processes affect the spectrum when the compactness parameter / is larger 
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than ~ 10 (Lightman 1982, Svennson 1987). When this condition is met pair cre- 
ation usually exceeds annihilation. Pairs lose most of their energy through Coulombic 
collisions. In some cases pair cascades, in which the relativistic pairs collide with the 
ubiquitous photons creating more pairs with less energy in the process and depleting 
the radiation field. If pair annihilation becomes important the 51 Ike V pair annihila- 
tion line may be present in the spectrum. The effects of pairs on the spectrum are 
a decrease of the (negative) spectral slope in the UV waveband due to loss of soft 
photons, and an increase in the hard X-ray - gamma - ray wavebands, due to the loss 
of pair producing hard photons. Pair production is traditionally calculated in the 
presence of relativistic electrons and UV photons of relative compactness //, for the 
electrons and l 3 for the photons. It should be noted that although pair production 
can happen without synchrotron or Compton processes, an input power-law spectrum 
is generally assumed for the electrons, mainly because most processes that produce 
enough highly relativistic electrons also result in a power-law photon spectrum. 

The atomic radiative free-free, bound-free, bound-bound and black body processes 
usually take place in ‘cold’(not completely ionized) plasmas. The features superim- 
posed on the continuum spectrum as a result are lines, edges and bumps in the IR 
through medium X-ray wavebands. The feature with highest photon energy is the 
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iron line-edge system near 7keV, which may have an absorption effect at up to lOkeV. 
(Absorption and emission features in the optical and UV regime are believed to be 
coming from regions which are far beyond the central engine. See Netzer 1990 and 
section 1.5.) 

Matter in local thermal equilibrium (LTE) may be present in the central engine of 
AGN, and may be responsible for the UV bump superimposed on the flat spectrum. 
This UV bump is thought to be black body emission from optically thick LTE plasmas 
with different temperatures. 

We shall limit our discussion to radio quiet AGN throughout the rest of the chap- 
ter. In radio quiet AGN most of the radiation is emitted in the EUV - X-ray wave- 
bands, and the X-ray emission is believed to originate in the central engine (see the 
following section). Radio loud quasars may emit significant amounts of radiation from 
their jets due to relativistic beaming effects (see appendix 1). As a result it is very 
hard to distinguish central engine emission from jet emission, x 

1.1.5 Variability of Emission from the Central Engines of 
AGN 

More than 50% of radio quiet AGN exhibit large amplitude flux variability in 
the X-ray waveband. Spectral variability is also present in many of these AGN 
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(Mushotzky, Done and Pounds 1993). Long term variability (years) is present in 
r> most sources, whereas at least half of the sources exhibit short term large amplitude 

variability (Green 1993, McHardy 1988). 

The flux variability (generally referred to as the light curve) is shown in figure 

O 

1-2 for three AGN. In the upper panel (l-2a) the light curve of NGC 4051 is shown. 
In the middle panel (l-2b) the light curve of NGC 5506. The light curve of MCG 
C -6-30-15, a Seyfert galaxy which is extensively discussed in chapter 6, is shown in the 

lower panel (l-2c). It is evident that, in the first and last cases, the total flux doubles 
on timescales of thousands of seconds. This change in flux is believed to originate in 

} 

the central engine. 
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Figure 1.2: Typical flux for two types of variability: a. XGC 4051 exhibits large amplitude 
'shot-noise’ time variability ; b. NGC 5506 exhibits an almost ‘chaotic’ time variability, c. 
MCG -6-30-15 has time variability which is in between a. and b. (see text). Most AGN 
have time variability which is similar to a. or c. (Adapted from McHardy 1988.) 
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The size of the emitting region is usually assumed to be smaller than R e = A tc, 
O where At is time in which the source is observed to double in flux. For incoherent 

radiation this condition is necessitated by causality. Events in a region of size R e 

t 

cannot be be correlated by perturbations that move faster than light. Since most 

r"\ 

of the energy is emitted in the variable part of the spectrum R e is the maximum 
possible size of the central engine. Since the shortest observed doubling timescale 
C goes to hours in Seyfert 1 galaxies (see table 1) the emission region can be assumed 

to be of order 10 13 cm in radius. The emission from this region is larger than ~ 10 42 erg 
sec -1 . As a result models that try to explain the light curve need to account for the 
enormous production of energy per volume, which is quantified by the compactness 
parameter in table 1. 

The Fourier transform of the light curves P = ^l/2x / 0 °° exp(iu)t) f (t)dt (usually 

referred to as the power spectrum) can be used to categorize the different types of 

variabilities. In the upper panel of figure 1-3 the power spectrum of NGC 4051 is 

shown, the power spectrum of NGC 5506 is shown in the middle panel and the power 

spectrum of MCG -6-30-15 is shown in the lower panel. Unlike many sources inside our 

Z> 

galaxy, there are no periodic components, which would appear as prominent spikes, in 
the power spectrum of AGN. The logarithm of the Fourier transform of a light curve 
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with random flux level (white noise) can be fit by a constant. The Fourier transform 
of flux from regions of random sizes and random flux levels is usually referred to as a 
‘random shot noise’. In this case the power spectrum log L(f) oc log(— 1/(/ 2 )), where 
/ is the frequency of variability, and the slope is —2. The flux of NGC 4051 is nearly 
a random shot noise. An almost chaotic intermediate logZ(/) oc log(— 1/(/)) power 
spectrum slope fits the Fourier transform of the light curve of NGC 5506. MCG 
-6-30-15 has an intermediate power spectrum slope of ~ —1.4. 
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Some models of the variability assume white noise which is filtered by matter in the 
line of sight. For example, when variable white noise flux is incident on matter with 
moderate optical depth the higher frequencies are filtered out. Spherical extended 
regions of size R e yield random shot noise from the original flat power spectrum 
(Sunyaev and Titarchuk 1980). This spectrum may be modified if pairs are present, as 
is expected in high luminosity to size ratio sources (Done and Fabian 1989, Lightman 
and Zdziarski 1987). Chaotic processes may also explain the observed power spectrum 
(Vio et al 1991), but no physical mechanism is yet known that may explain such 
processes. 

1.1.6 Geometrical Models of the Central Engines of AGN 

The most successful models of AGN include an accretion disk which channels 
accreted matter from the surrounding gas and stars onto a central super massive 
compact object in the central engine. In order to understand the observational success 
of this model one must look at models of the ‘outer regions’ (regions outside the central 
engine) of AGN. 

In figure 1-4 we see a simplistic, schematic description of the current geometrical 
model of the outer regions of radio quiet AGN. 
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Figure 1.4: This picture of AGN has an axial symmetry about the vertical center of the page. 
Radio quiet AGN are presumed to have molecular clouds in a torus configuration, electron 
scattering region, broad line region (the small open circles) and an inner accretion disk (the 
X shape in the middle). The observed emission therefore depends on the orientation of the 
observer with respect to the AGN axis. (Adapted from Mushotzky, Done and Fabian 1993.) 
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Our current picture of the outer regions of radio quiet AGN is: In the outermost 
regions lies the source of matter to be accreted onto the SMBH. This matter may 
originate in a rich stellar clusters, the inter stellar medium (ISM), or may be the result 
of cooling flows in clusters of galaxies (Blandford 1990, Rees 1990). We hereafter 
assume this matter to contain the usual cosmic abundance. Strong narrow lines are 
emitted by an unknown geometrical configuration of material with atomic number 
density 10 5 — 10 7 cm -3 , slightly more dense than the ISM, at distances of order 1 
Kpc from the central engine. This material is usually referred to as the narrow line 
region (NLR, see Netzer 1990). A little closer to the central engine molecular clouds 
form an axisymmetric ‘molecular torus’. The molecular torus may be responsible for 
the obscuring of parts of the more energetic parts of the spectrum, and is probably 
responsible for most IR emission. Closer to the center is the Broad Line Region 
(BLR), the geometry of which is uncertain. The broad line region clouds are the 
open circles near the center of figure 1-4. From the observed width of the optical 
lines we know that the BLR matter is moving fast. Other observations reveal the 
filling factor is small, and the number density of matter is 10 8 - 10 lo cm” 3 in the BLR 
(Netzer 1990). Evidences for a super massive compact object and a superluminous 
central engine are mostly found from observations of NLR and BLR emission lines 
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(Netzer 1990). 

Closer to the central engine a geometrically thin accretion disk, which is optically 
thick, is usually considered to be present. In figure 1-4 the disk is represented by 
the central dark bow-tie shaped region. Matter in this disk is accreted onto the 
SMBH. This disk is generally considered to emit a significant portion of the total 
luminosity emitted by AGN in the UV band (the so called ‘UV bump’), and may 
also be responsible for the ionization of clouds in the BLR. The disk configuration is 
most efficient in the outwards transfer of angular momentum, which allows for high 
accretion rates. Angular momentum is viscously transferred from accreted matter to 
outflowing particles and radiation. 

Since the model for the outer regions cannot explain the X-ray emission and 
variability a different, more elaborate, model is needed for the central engine. The 
model of the central engine must take the accretion disk into account. 

Calculations show that the thin accretion disk is thermally unstable near the 
central engine because the heating due to viscosity exceeds the radiative cooling (see 
Lightman and Eardley 1974, Thorne and Price 1975). As a result, a geometrically 
thick disk which can be optically thick (if supported by radiation pressure) or optically 
thin (if supported by ion pressure in a two-temperature plasma, see Shapiro, Lightman 
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and Eardley 1976; Rees et al. 1982) may form. The central engines of AGN are 
considered to be the inner part of such disks. The electron temperatures of these disks, 
which may exceed 10 9 K, and the magnetic fields which are of order 10 3 gauss, can 
account for the synchrotron and Compton processes that can produce the observed 
spectrum. When pair-cooling is important the central regions of the geometrically 
thick disk may collapse into an assembly of clouds, an annulus, or thin disk. The effect 
of pair plasmas on optically thin and optically thick disks is discussed in several papers 
(Begelman Sikura and Rees 1987, Tritz and Tsuruta 1989, Kusunose and Takahara 
1989). 

In figure 1-5 we show the schematic description of the central engine, which as- 
sumes a disk configuration. Spectra for the disk corona model were calculated by 
Lightman and White 1989, George and Fabian 1990, Tritz 1990, Tsuruta and Kellen 
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Figure 1.5: The disk model for the central engine of AGN includes an outer cold disk which 
flares up to form a thick torus of very hot plasmas (10 9 -10 12 K) in its inner regions, and 
collapses back to a thin disk configuration near the SMBH. 
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An alternative model of the central engine is the ‘disk corona model (Svennson 
r and Zdziarski 1994). In this model the disk is not disrupted because of pressure from 

the corona in which most of the radiation is produced. The radiation in this model 

is mainly through Comptonization of the UV disk emission. 

O 

Many models attempt to reproduce the various AGN light curves while assuming 
the existence of an accretion disk in the central engine. One model which suggests 
{$ that the magnetic field reconnection at the surface of an accretion disk may result 

in an acceptable power spectrum (Pudritz and Fahlman 1982, De Vries and Kuipers 
1989, see also remarks in chapter §7). Another disk model include randomly appearing 

o 

uniformly distributed hot spots, which are related to magnetic flux tubes and may 
appear on the surface of Keplerian accretion disks with different apparent intensities 
due to the Doppler effect (Abramowicz 1991). 

An alternative model for the central engine which does not include an accretion 
disk is presented in section 3 of this chapter. 

1.2 INTRODUCTION TO X-RAY OBSERVA- 
TIONS OF AGN 

1.2.1 Why X-rays? 

The central engines of radio quiet AGN are most readily observed in the X-ray 
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waveband. Radio quiet AGN are observed throughout the electromagnetic spectrum, 
but most of their energy is emitted in the UV, X-ray and 7-ray wavebands. Recent ob- 
servations by the OSSE and COMPTEL instruments on board the Compton gamma 
ray observation satellite have shown a high energy end cutoff in the hard X-ray - soft 
7 wavebands of most radio quiet AGN, above a few hundred keV (Ghiselini et al. 
1994 and references therein). The IR observations of the same sources reveal a cutoff 
at 10~ 3 eV. The dominant parts of the spectrum are therefore in the UV - soft 7-ray 
wavebands. Unluckily the EUV part of the spectrum is heavily obscured by ISM 
in our galaxy, thus making observations extremely difficult. The X-ray waveband is 
therefore the most dominant part of the spectrum that can be observed with present 
technology. 

Evidence that most of the EUV and X-ray emission originate in the central engine 
mount. Studies of the BLR show that an internal ’hard’ (photon energy higher than 
few tens eV) source is irradiating the BLR matter. It is generally accepted that in 
radio quiet AGN there is no strong ‘contamination’ of the X-ray spectrum from the 
central engine with emission from strong relativistic jets (although faint small scale 
jets are observed in Seyfert galaxies, Wilson 1985, Mushotzky, Done and pounds 
1993). The IR spectrum of these sources, which is generally thought to originate 
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in dust from the molecular torus, is not variable (Done et al. 1991). The optical 

p continuum, which may originate in the thin accretion disk, cannot dominate because 

it covers less than a decade of frequencies (see figure 1-1). 

The best evidence that most of the X-ray flux originate in the central engine 

C 

comes from an argument about large amplitude flux variability. Timing studies show 
that the size of the X-ray emitting region is much smaller than the BLR region (see 
C sections 1.4, 1.5). Whereas in other wavebands the variability of radio quiet AGN, if 

existing, is on timescales of years and with very small amplitude (see Mushotzky Done 
and Pounds 1993, for example), the large amplitude variability of X-ray emission is 

r. 

rapid (doubling the flux in hours, see for example section 1-4 figure 1-2). X-rays 
are, therefore, the window through which we are able to collect most data about the 
central engine. 

In the observational part of this thesis we concentrate our attention on Seyfert Is 
nuclei, a large subgroup of radio quiet AGN. Seyfert Is show significant variability 
amplitude over a short time scale in the X-ray waveband and may therefore provide us 
with more information than QSO. Since the time scale of variability of QSO is years 
long term observations from X-ray satellites are not practical. Seyfert 2 galaxies are 
excluded because there is evidence of large absorption columns in the line of sight 
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(probably due to the molecular torus, see Antonucci 1993). These large absorption 
columns are probably responsible for the reprocessing of the original central engine 
radiation. 

1.2.2 Large Area Counters, X-ray Telescopes and Other 
Technologies 

X-rays cannot be observed through Earth’s atmosphere. Balloons, rockets and 
satellites are therefore used in order to go above the atmosphere. Observing from 
very high altitudes may not be advantageous for several reasons. One reason is the 
increased cost of the observation. Another is that Earth’s upper atmosphere actually 
operates as a shield against energetic cosmic rays which may excite nuclei and high 
Z atoms in the satellite, resulting in internal background due to atomic and nuclear 
decay. Most X-ray missions are therefore launched into the so called low-Earth orbit, 
a few hundred kilometers above Earth’s surface. 

Very few X-ray photons from extra galactic sources are collected by the detectors 
on board X-ray satellites. The internal background may therefore be much larger 
than the signal, requiring a sophisticated array of photon rejection techniques. The 
X-ray background radiation, which is poorly understood, should also be rejected when 
pointed observations are performed. Large area counters, which collect the largest 


'J 



number of photons, are the most sensitive to these types of noise. 

In X-ray telescopes, on the other hand, the detector volume is small, and the 
internal background is therefore subdued. Photon statistics, on the other hand, are 
very poor in these detectors, and the sensitivity to flux variability is limited. 

1.3 AGN IN THIS THESIS 

1.3.1 Theoretical Treatment of AGN in This Thesis 

In the first part of the thesis we explicitly calculate the effects on the spectrum of 
the ‘two component cloud model’ of the central engine, which under certain conditions 
is a viable alternative to the disk and disk corona models. The two component cloud 
model was first suggested by Guilbert and Rees (1988). Guilbert and Rees and Celotti 
Fabian and Rees (1992) noted that moderate compactness values may result in high 
optical depth r = noR , where n is the number density of the accreted matter, a 
is its cross section for absorption and R is the size of the central engine. Density 
perturbations thus grow as a result of radiation or magnetic pressures. When the 
density is increased cooling in the dense regions is increased. The cooler matter has 
even higher optical depth, and the added radiative (or magnetic) pressure increases 
the density even further. The process terminates when the denser regions become 



thermalized and self pressurized. These ‘cold’ denser regions, which may be in the 
form of plasma sheets and filaments, are referred to as ‘clouds’ in this thesis. Typical 
temperatures are 10 9 - 10 12 K for the hot component and 10 5 - 10 7 K for the cold 
component clouds. It is not at all clear that these clouds can form an accretion disk. 
Since the clouds are embedded in a hot plasma this model is referred to as the two 
component cloud model, or ‘the cloud model , in this thesis. 

In figure 1-6 the general picture of the cloud model is presented. 
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Figure 1.6: The cloud model (or: ‘two component cloud model’) for the central engine may 
include an outer cold disk which flares up to form a thick torus with two phases: very hot 
plasmas (10 9 -10 U K) in one phase, and a colder dense plasmas (10 -10 K) in the other phase. 
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The clouds may manifest themselves in emission, absorption and reprocessing 
spectral features which characterize cold matter. Recently such features were discov- 
ered in X-ray observations of AGN (e.g. Mushotzky, Done and Pounds 1993). These 
clouds may also be responsible for the UV bump which is frequently observed in the 
spectrum. 

In chapter 2 we present an equation of state for the cloud model. The radiative 
transfer equations for clouds in chapter 3, which utilize the equation of state, results 
in spectra which are shown to generally fit the observed spectra of radio quiet AGN. 
In chapter 3 we also describe the conditions under which the two component cloud 
model is viable. 

In the second part of the thesis we present an alternative model for the variability 
of the central engine. A problem with the disk and disk-corona models is the a 
posteriori nature of the explanation of variability. We develop an alternative, self 
consistent, model for the variability of the central engine in chapter 4. We investigate 
the dynamical properties of the central engine, and we also consider the possibility 
of shocks. We show that shocks, which readily form in the central engine, may be 
responsible for the observed light curve and spectral variability. Shocks are also 
considered to be efficient in transferring angular momentum, accelerating particles 
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and the creation of the cold and hot components of the cloud model. The importance 
of shocks in the calculations of the stability of optically thin disks is discussed. 



1.3.2 X-ray Observations in this Thesis 


In the third part of the thesis we present results of the analysis of X-ray data from 
several one day observations of the Seyfert 1 galaxy MCG -6-30-15 over a period of 
five years. This extensive analysis was designed to reveal the possible structure of the 
central engine in this source. These results are particularly interesting in the context 
of the theoretical part of this thesis. 

The observational part of our work is discussed in chapters §5 and §6. In chapter 
§5 we describe physics of proportional counters which are used in the detection of X- 
rays. The detectors on board the Ginga and ROSAT satellites, with which the data 
analyzed in this thesis was collected, are also described. In the last part of chapter §5 
we describe the methods by which background is subtracted and the method by which 
the spectra are selected and sorted. The analysis of data from the Seyfert 1 galaxy 
MCG -6-30-15 is described in chapter §6. The analysis of spectral and flux variability 
and the possible explanations in term of simplified models are also presented in that 
chapter. Implications of the observations, and other X-ray observations of Seyferts, 
are discussed in chapter §7, the last chapter. 



Chapter 2 

EQUATION OF STATE FOR 
THE CLOUD MODEL 


In this chapter the equations of state for clouds in the central engines of AGN are 
presented. In the first section we present the theoretical and observational arguments 
for a two component configuration. In the second section we develop the equations 
for a simplified spherically symmetric case. In the third section we solve the equations 

V 

for a stationary case in which the physical parameters are distance dependent, and 
compare the results with other solutions for the central engines. 

2.1 REASONS FOR THE TWO COMPONENT 
CLOUD MODEL 

2.1.1 Theoretical Reasons for the Cloud Model 


The approximate radiation density and matter density in the central engines of 
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AGN can be found from observations. For radio-quiet AGN the observations indicate 
that the densities, and the resultant compactness (see table 1), may induce strong 
coupling of matter and radiation. This coupling results in rapid local cooling, and 
the creation of two components: the radiative ‘hot’ component (ksT ~ me 2 ) is 
dispersed throughout the central engine, and the denser ‘cold’ (not completely ionized) 
component, which reprocesses the radiation, form clouds. 

The radiation and matter densities, pressure and temperature of the central en- 
gines of AGN can be determined by using three observed quantities: the luminosity, 
the doubling timescale and the absorption column, which is the number per cm 2 of 
atoms in the line of sight to the inner parts of the central engine. 

We shall limit our discussion to a large sub-group of AGN, radio quiet quasars 
(QSO) and Seyfert nuclei, in which the central engine is most readily observed. 
Seyferts and QSO may be the ‘cleanest’ example of AGN as far as exploring the 
central engine is concerned, and are also the most abundant (Woltjer 1990, appendix 
I). Another advantage of limiting consideration to this subgroup is that some QSO, 
which were identified as point sources when first observed, have already been resolved 
and found to be galaxies which harbor an active center. A large range of observational 
parameters is therefore available for comparison with our model. 
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The AGN flux from QSO and Seyfert 1 nuclei is observed to double (in X-rays) in 
hours (some Seyferts) or days (quasars), which implies that their size is much smaller 
than that of the host galaxy. If the variability amplitude is large, the size of the 
source is determined by the doubling time scale 

Rx=cAt , (2.1) 

where Rx is the size of the X-ray source, c is the speed of light and At is the observed 
shortest time scale of large amplitude variability. In AGN this time scale is usually 
taken to be the minimum interval in which the source doubles its output. For the 
observed variability timescales the size which is found to be Rx = 10 12-15 cm is of 
order of the size of the solar system at the large end. The size of the SMBH must be 
smaller than the size of the X-ray source but cannot be much smaller, or else matter 
will not be efficiently captured and accreted. Assuming that a SMBH is responsible 
for the accretion of the emitting matter and that the emission comes from ~ 10 , 
where R g = GM/c 2 , and the mass of the SMBH is found to be of order 10 6_9 A/q, 
where Mq = 2 x lO^gr is the mass of the sun (Blandford 1990). The efficiency of 
accretion onto black holes varies from ~ 0.06 for accretion onto a Schwarzschild black 
hole to ~ 0.4 when the accretion is onto a fast rotating Kerr black hole. Most of the 
energy is released near the horizon of the black hole, at ~ 10f? a (Blandford 1990). 
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Radio and optical observations have found the central source to be smaller than 
~ O.lpc, .which, however, is still more than three orders of magnitude larger than 
the suggested X-ray region. This is an upper limit on the size. The smaller sizes 
are also inferred from the mechanisms for the production of jets, where these are 

0 

observed in other types of AGN (using certain unification of schemes of AGN). Future 
observations by missions such as the COSTAR on the HST and the advanced X-ray 
facility (AXAF) will determine weather the smaller size is valid for these sources. 

Observed values for the total luminosity vary from 10 41-44 ergs sec -1 for Seyferts 
to 10 44-47 ergs sec -1 for QSO, requiring an extremely efficient radiation mechanism. 
There is an inherent uncertainty in the above numbers, due to lack of observations in 
EUV regime in which the emission is generally expected to be the strongest. However, 
recent observations by the European ROSAT satellite in the soft X-ray waveband, 
and by the international UV explorer (IUE) in the near UV band, have confirmed 
the existence of excess of radiation in the EUV band in many cases (see also section 
§ 2 . 1 . 2 ). 

The supply of matter to the central engine must exceed 

I 

M = - j, (2.2) 

7]C l 

where 77 is the fraction of accreting matter turned into energy and L is the observed 
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luminosity. The rate of accretion is therefore M ~ (10 2O— 23 /r/)gr sec -1 for Seyferts. 

A primary reason for a two component medium is the physical availability of pres- 
sure that enhances perturbations. For example, when the source radiation pressure 
exceeds the gravitational pull the overly dense regions are compressed by radiation 
pressure. These compressed regions emit synchrotron and Compton radiation and 
may cool. The cooler regions have higher cross-section for Compton scatterings, and 
are therefore compressed and cooled even more. When the matter is sufficiently cool 
the bremsstrahlung and recombination radiation, which depend on the square of the 
gas density, dominate. Further cooling and density increase stop when the matter 
reaches local thermal equilibrium (Guilbert and Rees 1988, Ferland and Rees 1988). 
Although radiation pressure is easily understood, the main reason for this local den- 
sity increase may be twisted magnetic fields, probably in a configuration similar to 
that observed in quiescent solar prominences (such as in Priest, Hood and Azner 
1989), which channel cold matter into low field high gradient regions. 

Observations imply that cold matter is close to the center, or within the central 
engine. The above theoretical argument for the existence of cold plasmas in the 
central engine comes from measuring source compactness. Sources which are near 
their Eddington limit are optically thick if the accretion is locally homogeneous and 



isotropic. When the Eddington luminosity is reached matter is accelerated outward, 
because the radiative force on each particle equals the gravitational force 

LEddCT _ GMprnj^ ^ ^ 

47rR 2 C R 2 ’ K ' 

where L E dd is the Eddington luminosity, <Jj is the Thompson cross section and y.m H 
is the mass of an average cosmic abundance atom. Accretion may continue if it is not 
isotropic, depending on the configuration. 

This result may indicate that a thin accretion disk is the preferable geometrical 
configuration for the central engine. However, it has long been known that the spec- 
trum and light curve of radio quiet AGN cannot be explained by the standard thin 
disk model alone, and that an optically thin hot relativistic plasma in the vicinity of 
the compact object may be responsible for a large portion of the continuum X-ray 
emission (Price k Thorne 1975, Shapiro Lightman k Eardley 1976). This region may 
form when the inner regions of a thin disk inflate due to thermal instability (Lightman 
k Eardley 1974). Models with thin cold disk and hot corona are discussed in chapter 
7. We show that high compactness sources must have an optically thick component 
somewhere in the hot relativistic plasma (Guilbert and Rees 1988, Lightman and 
White 1988, Sivron and Tsuruta 1993). The location in which the cold component 
forms is arbitrary, and it is therefore imperative that we shall start our study of the 
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central engine with a general form of cold matter distribution, not only an a priori 
disk configuration. 

2.1.2 Observational Reasons for the Cloud Model 

»-x There are three types of observations that imply the existence of cold matter 

in the central engines of AGN: The spectrum of AGN continuum radiation may be 

described as a combination of of hot and cold components spectra; The temporal 

variability may be described as a typical hot component variability which is filtered 

by the cold component; Multi-frequency observations of short term variability may 

provide evidence for the existence of cold matter inside the central engine. 

The spectrum of a typical radio quiet AGN is almost ‘flat’ (see figure 1-1), which 

is typical of hot, optically thin, matter. This property was the reason for the variety 

of theories that explained the spectrum as one power-law from one region. One 

such theory is the synchrotron self-Compton model (SSC, see Blandford 1990), in 

which a power-law spectrum from the synchrotron self absorption edge at v ~ 10 12 

sec -1 to the pair-plasma ‘break’ at v ~ 10 2 ° sec -1 (see Ghisellini and Haardt 1994 

and references therein) gives rise to features such as a possible hard x-ray tail, by 
J 

self Comptonization of the same synchrotron radiation. In this model the source 
of emission is an almost optically thin relativistic plasma of temperature 10 K, 


which is located in, or near, the central engines of AGN. In another mechanism, 
usually referred to as the Compton mechanism, pairs and electrons in a hot plasma 
component in the central engine scatter photons from a ‘cold’ (not completely ionized) 
component, which results in the creation of a power law spectrum. 

Recent observations in the UV and soft x-ray parts of the spectrum have revealed 
that a large portion of the luminosity in these wavebands is better characterized as 
thermal spectrum, typical of a cold, optically thick, matter. The overall portion of 
the total luminosity that can be characterized as thermal UV bump is: 

0.3 < - < 3, (2.4) 

Lnt 

where Lx is the thermal luminosity and Lnt is the power-law non-thermal luminosity 
(Lightman and White 1988). 

Another reason for considering a two components cloud model comes from studies 
of the random shot noise and chaotic variability which is typically observed in AGN 
(see section 1-5 in chapter 1). The power spectrum can be best fit with a 1// slope 
in some strong AGN, or with 1 / f 2 for some very weak AGN and galactic black-hole 
candidates (McHardy 1988; Green 1993, and references therein). The variability in 
high luminosity sources can be accounted for if there is an increase in cold matter 
which reprocesses the original random shot noise in these sources. The slope of the 
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power spectrum can also be flattened because of pairs, which are more copious if the 
(\ luminosity to size ratio is large, and may increase the rate of cold component matter 

production (see also chapter 4). 

Multi-frequency observations have not yet revealed the existence of cold matter 

O 

in the central engine. In the recent multi-frequency observation campaign, organized 
by Prof. S. Tsuruta of Montana State University, it was found that at least in one 
C example less than 5% of the luminosity in the IR and optical parts of the spectrum 

vary simultaneously with the mostly power-law X-ray emission (Done et al 1991). 
This observation may rule out the SSC model, or imply the filtering out by clouds of 
the optical - IR components of the spectrum (Celotti, Fabian and Rees 1992). Future 
observations with better statistics and improved spectral techniques may reveal better 
evidence for cold matter in the central engine. 

2.2 THE EQUATION OF STATE 

In deriving the equations of this section the general approaches of Guilbert and 
Rees 1988, Lightman and White 1988 and Celotti Fabian and Rees 1992 are followed. 
The observables used in our framework are the X-ray luminosity Lx, the mass of 
the SMBH M and the total number of cold absorbing atoms in a cm 2 column in 



the line of sight N H (the absorption column density). In the computer code we used 
the fraction of the total Ax solid angle available for light to freely escape the central 
engine, 0 <C ^(AT/f, A/, £x) ^ 1 (the covering factor) as a variable. In this chapter 
our attention is confined to accretion-powered AGN radiating close to the Eddington 
limit. 

We first describe the isotropically accreting matter. The total luminosity in equa- 
tion (2.2) determines the isotropical accretion rate and local density 

M AxR 2 Tt a \jg YYl fj t\j & (2.5) 

where v acc is the velocity at which matter is accreted onto the SMBH. The accretion 
rate can be also written as 

M = Ax R 2 n avg mf{ (v acc /vjf) \j2GMjR (2.6) 

where n avg (R) ~ Pavg/mH is the average number density of the infalling matter. 
Since the cosmic abundance presumed is mainly made of hydrogen with mass m# we 
assume this to be the mass of all infalling particles. (The abundance can easily be 
changed by replacing Mh with an arbitrary pm H - We choose p - 1 to be consistent 
with Guilbert and Rees 1988 and Celotti Fabian and Rees 1992.) The free fall velocity 
vjj is assumed to be Newtonian. As a first order approximation the accretion rate, 


which is the solution of the spherically symmetric continuity equation, is assumed 
to be independent of R, the distance from the SMBH. This approach renders our 
treatment stationary. Dynamical considerations are discussed in chapter 4. 

The minimal optical depth of isotropically infalling matter is determined by 




(2.7) 


where k e , is the opacity for electron scattering. Substituting equation (2.6) in equation 
(2.7) and integrating over the central engine we get: 


r > 0.4 


R L / 
4 ttR 2 t}<? \ 


vjf\ nr 

v acc J V 


■ ( 2 . 8 ) 


where R here is the size of the central engine. Equation (2.8) can be written in a 


more convenient form, 


^ ,/! fe) - (2 - 9) 

where t acc , the accretion timescale, is assumed to be proportional to £//, the free 
fall timescale. This dependence shows that geometrically extended central engines, 
such as those expected in the case of extended two temperature ion torus (Rees et 
al. 1982), may survive without transition to two components. The optical depth r 


would exceed 1 near the center when the source is radiating close to L e< M (Guilbert 
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and Rees 1988). Thus the original power-law emission, which is probably emitted 
from the inner 10i? 9 (Blandford 1990), will inevitably be shrouded and reprocessed 
by the infalling matter. Consequently the final spectrum deviates significantly from 
a simple power-law even if the luminosity is an order of magnitude smaller than the 
Eddington luminosity (Guilbert and Rees 1988, Lightman and White 1988). The 
slight inverse dependence on the distance is the result of the crude approximation 
adapted, and is especially inaccurate for large optical depths which include processes 
other than electron scattering. 

Using the same arguments as the ones in equation (2.7) (where we use (2.15) for 
the density) and integrating radially over one cloud, the optical depth of a cloud in 
the two component configuration is calculated. Each of the cold plasma clouds has 


optical depth 


[ Vmax , (k u + kkn ) 3 

Tlcld = / dv- 

JVr, 


(-L _ 


( 2 . 10 ) 


/ ^mir 


in 167r 2 77n c /d/^ d (u acc /u//)\/2GM \y/R y/R + hid) 

where k K n and the Klein Nishina and bound-free opacity, respectively, replace 
the electron scattering opacity. The radius of one cloud is n c (d is the number 
density of clouds at the particular radius, and v is the frequency. From radiative 
transfer (see chapter 3) one can see that the total high energy flux coming out of this 
region is roughly L(out) = Lexp^—T\ c i,iRHcid\ A more accurate optical depth, such 




as the one presented by Ferland and Rees 1989 and Netzer 1993, includes the atomic 
properties and the Einstein Coefficients. 

Because of fast Compton and faster bremsstrahlung and bound free cooling most 
of the shrouding infalling material should be relatively cold. The minimum temper- 
ature can be derived for thermal matter which emits as much energy as it radiates: 
M 2 dd a B T* = (1 - A)L( 1 - exp(-T lc{d ))exp(-f(R))/(hrR 2 ) from which we get 

Tdd = ((‘ - A) (' “ e~ /(fi> ) 1/4 , (2.11) 

where A is the flux averaged albedo (per cloud) and f(R) is the covering factor at R 
(see equation 3-12 for f(R) and 3-10 for the albedo). An equality is reached if n ( A, the 
cold cloud’s number density, is high enough so that the cooling time hound- free « 

hremj ^synchro— Compton fpatrs, where Abound— free 1® the timescale for Cooling due 

to matter that is not completely ionized. 

Such densities can result from a confining pressure supplied by an equipartition 
magnetic fields, P mag netic = B 2 /16n ~ GM/R , as in Rees 1987, (or, in a disk con- 
figuration, from the differential gravitational pull of the black hole). A two phase 
structure with dense cold clouds and hot relativistic component forms when the cool- 
ing time in the dense phase is shorter than the infall time scale by a factor which 
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accounts for the ratio of photon energy to particle energy 

h rem ,(T) < (GM/^y^ikT^/rjm^ 2 ), (2.12) 

where hrems = oJ 1 (kT c id/m e c 2 ) 1 ^ 2 l/(n e (TTc). This cooling argument may also be 
explained as follows: since in a typical infall time scale a particle in the hot phase 
gains energy of order GMm p /R and since at the same time the cold phase is cooling 
via bremsstrahlung in which every collision releases a photon of energy which exceeds 
kTgg, there must be at least enough such photons released such that no energy is 
gained. The above relation will result in a very high density. 

By equating the thermal pressure with the confining magnetic virial pressure we 
get (Guilbert and Rees 1988, Lightman and White 1988): 



Note that the inverse of this ratio is approximately the filling factor. Also note 
that the radiation pressure is proportional to the virial magnetic pressure, but that 
dynamical motions are than required. (Equation (2.13) presumes the gas equation for 
the hot component. This assumption is valid for the stationary equations. In cases 
where the total energy of the hot component doubles in seconds this assumption may 
not be valid. See also chapter 4.) 


j 



o 


53 


Combining equation (2.12) with (2.13) , and using (2.4) and (2.6) for the average 
density and the Eddington luminosity, we get that when the inequality 

(^y'^Ly' (R\ (u 


OTi 


(£)" £) (t 


-2 


< 1 


(2.14) 


is fulfilled a dense phase will stay cold. It is evident that this inequality holds for 
near Eddington luminosities. 

Using the same arguments that lead to equation (2.5), the isotropic accretion of 
cold component clouds should obey: 


dMcia _ 47T fl 2 Tim 4 irQj n ih m H y Ml 

dt 3 

where v c m acc is the accretion velocity of clouds. 


(2.15) 


If the requirement that the clouds hold a significant fraction of the infalling mass 
is fulfilled a condition on the size of the clouds can be found. We rewrite equation 


(2.16) as 


dM L£ 4 ttR 2 n c i d 4t l* ld n th m H v a 


(2.16) 


dt rjc 2 3 

where £ is the fraction of the mass in the clouds, and M = const(R). This condition 

will be fulfilled if a significant fraction of M in the form of clouds does not fall into 
the equatorial plane to form a geometrically thin accretion disk, when the clouds are 
small and dense enough so that the time scale for cold matter to fall onto a disk 
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is longer than the typical infall time scale Rjvjj. We take the typical collision time 
tcoii ^ l/{ n cidftlcid v cid) to be the time scale of falling into a disk. From this requirement 
we obtain the following relation between n c id and Idd' 

n ddldd < (kRY 1 - (2.17) 

The cloud velocity v c id is uncertain, but it may be assumed to be Keplerian. Plugging 
equation (2.16) in (2.17) and using the definition of compactness in section 1.1 of 
chapter 1 one get 

nthldd > 10 24 cm 1 - (2.18) 

We see that semi spherical accretion in radio quiet AGN may take place for clouds of 
sizes less than I c \a ~ 10 8 cm and atomic number density n c /<* ~ 10 16 which is expected 
in the central engine. A typical cloud of size 10 6 — 10 8 cm is much smaller than the 
size of the central engine in radio quiet AGN. The above consideration is a simplified 
form of the macroscopic equation of motion for clouds which is further considered in 
chapter 4. 

When the condition of equation (2.18) is not met the configuration is a disk full of 
clouds, or clouds which collapse to form a flat thin disk with sheets and filaments that 
still hover above and below the disk. In such case equation (2.16) includes another 


j 



parameter for the fraction of accreting matter in the clouds. 


2.3 RESULTS 

By simultaneously solving our basic equations, (2.10), (2.11), (2.14) and (2.17) 
the radial distribution of our characteristic parameters is obtained. For the spherical 
accretion of inhomogeneities we assume functional R dependence for the following 
quantities: 


Tlth — ^th(,Rmax) 

( R f 

(2.19) 

T c ld = Tcld(Rmax) 

( r r 

'Rmax ' 

(2.20) 

T^dd = W'ddi'Rmax) 

& 

(2.21) 

^cld ~ ( 

f r y 

^Rmax ' 

(2.22) 


Radial dependence is justified for high enough accretion rates, such that r > 1 in 
equation (2.7), and small enough clouds, such that (2.18) is fulfilled. The maximal 
radius, R ma x, is the outer boundary in which the cloud model works. The reason that 
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the clouds, sheets and filaments have distance dependence physical properties, is that 
magnetic stresses, or the shocks, which probably result in the formation of the cold 
phase, are distance dependent (see also chapter §6). Equations (2.19) - (2.22) give us 
four unknown exponents. 

These exponents can be found if the covering factor and albedo are small and the 
optical depth per cloud is large. The radial dependence of the temperature assumes 
small covering factor, and thus it is easily seen from equation (2.11) that T c id oc R~ x t 2 , 
or A = 1/2. Furthermore, if the virial dependence of the confining pressure holds, as 
in equation (2.13), we find that murith oc f? _5 ^ 2+A , which gives an inverse square R 
dependence for the number density of cold plasma in the clouds with small covering 
factor. This result is encouraging. In a typical Seyfert galaxy the density needed for 
LTE in the central engine is approximately n t h = 10 16 cro -3 for a Bolometric luminos- 
ity of 10 A4 erg j sec, which is the minimum needed for there to be an effective clouds 
coverage. One may compare this density to the broad line region clouds at ~ 10 4 R g , 
which have a density of n t h ~ 10 lo cm -3 in Seyferts. From this observational result 
we also get n t h oc R~ 2 . When we include the theoretical results of Celotti, Fabian 
and Rees 1992 which require that the number density in a cloud near the SMBH is 
approximately n t h = 10 16-l8 cm -3 for an absorption column of Nh < 10 21 cm -2 , such 
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that the IR-optical wave bands are blocked we find that n t h oc R~ 2 — R~ 2GG . We like 
to stress, however, that this does not imply that the BLR clouds and the clouds in 
the central engine a part of the same physical mechanism. Using equation (2.18) the 
dependence of the clouds’ size exponent have i ~ 1 — 2 and the cloud number density 
has the exponent x = 3 — 4. 

With the above results we are confident that the equations above can be used in 


the radiative transfer code. 



Chapter 3 

THE RADIATIVE TRANSFER 
FOR THE CLOUD MODEL 


The radiative transfer equation for a two component cloud model is presented 
and used for the calculation of output spectra in several wavebands. We assume 
that these clouds have many similar properties, and develop a macroscopic radiative 
transfer equation. The radiation transfer problem is written for macroscopic bodies, 
such as a multitude of clouds, and an optically thin hot component. The spatial 
distribution and physical properties of clouds roughly follow the equations of chapter 
2. Possible solutions of the radiation transfer problem are outlined. 
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3.1 RADIATIVE TRANSFER FOR THE TWO 
COMPONENT MEDIUM 

In this section the radiation transfer equation is described. We compare absorption 
and emission by clouds to the well known atomic absorption and emission processes. 
We than combine absorption and emission to form the full radiative transfer equation. 

3.1.1 Absorption and Emission by Macroscopic Inhomo- 
geneities 

Absorption by atomic processes and absorption by clouds are similar. The ab- 
sorption due to atoms with cross section cr v is dl„ = — a v nlds , where / is the specific 
intensity (Rybicki and Lightman 1979), the absorption coefficient is a v = cr v n, and 
the solution for constant a is / = I(s = 0) e -a *. The absorption by cold clouds can 
be expressed as 

dI Vl n = (-M(fl) n dd I v ds, - (3.1) 

where A(O) is the effective area of a cloud and the subscript cld is used indiscrim- 
inately in this chapter for clouds of arbitrary and specific shapes. The area of an 
opaque spherical cloud, Am = tcRm, is the same as its cross section. (Optically thin 
clouds have frequency dependent smaller cross section.) For spherical opaque clouds 
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we therefore have a typical extinction of the intensity: 




(3.2) 


In figure 3-1 we see the expected distribution of absorbing clouds in AGN and the 


resultant outcoming radiation. 



Source of Power law 



Outcoming radiation 


Figure 3.1: Absorption by optically thick clouds 
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Atomic materials and clouds both have direction dependent, and direction inde- 
pendent types of emissions. The emission generated by atomic material is 


h / = 


dE 


dVdVldudt' 


(3.3) 


The intensity added to a beam traveling a distance ds through a volume element of 
unit area in an inhomogeneous material is 


dl u = jvn cli ds, (3.4) 

where j v is the emission coefficient for one cloud, and an exponential gain from the 
medium is possible, if j u oc /„ (see, for example, appendix 1). Atomic emission may 
be of the type 

. _ e 2 U; 4 A 2 COS 2 # __ PLarmor 

^ 3c 2 47r 4 tt 

where A is the amplitude of oscillations of a dipole. Thermal emission is an example 
of an isotropic emission. A black-body (KirchofF’s - law) emission is: 

2hu 3 /c 2 

j u n<7v t xp^hi/jkT) — l’ 

which is the result of equilibrium of absorbed and emitted radiation. 



(3.5) 



The emission by clouds, which is not necessarily isotropic, has the form: 

^ nc,d = {dv dn It du) ’ ^ 

where j„ is the emission of one cloud, and is the number of clouds per unit volume. 
(This number will usually be very small, of order 10 -3O cm -3 ). We define the emission 
per cloud to be j u , because of the simplicity in adding cloud emission spectra to the 
computer code. For reflecting clouds j v is direction dependent. In this thesis, for the 
case of optically thick matter, we simplified this problem by using the ‘moon phases 
method' (see appendix 2). For optically thin matter a Monte-Carlo simulation for one 
cloud will give an angle-dependent emission (see, for example, Bond and Matsuoka 
1994 and George and Nandra 1994). 

In figure 3-2 we see the expected emission from an optically thick cloud in the 


central engine of an AGN. 
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Figure 3.2: Compton reflection and thermal emission from optically thick clouds 
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3.1.2 Radiative Transfer in a Two Component Medium 


Several simplifying assumptions have to be made in order to write a macroscopic 
radiative transfer equation. The radiative transfer equation for clouds is macroscopic 
in the sense that the absorbing and emitting elements are macroscopic bodies. In 
our code the non-thermal source is represented by a point source that emits a power- 
law spectrum from the innermost region of the central engine. This assumption is 
justified by the results of section 1.2 in chapter 2 and section 2 of chapter 4. The 
model also has to satisfy the variability requirement imposed by the absence of IR- 
optical variations during X-ray flares from Seyfert nucleus NGC 4051 (Done et ai. 
1991). Therefore, we assume that the IR-optical portion of the primary radiation 
is effectively absorbed by the extremely small ‘cloudlets’ within 10/2 3 , which are 
transparent to the higher energy photons. Celotti Fabian and Rees (1992) predict 
the existence of such cloudlets, which we shall call ‘the CFR cloudlets . In our cloud 
model this modified spectrum is used as the input radiation, to be reflected and 
reprocessed by an assembly of cold, optically thick clouds. (The energy from the IR- 
optical wavebands is added to the UV band.) The inner boundary in our calculations 
is therefore ~ 10 R g . (We note, however, that in models with pair plasmas the intensity 
of optical-IR radiation from the central engine is inherently small, and CFR cloudlets 


are not necessary.) 

The clouds in this thesis are distributed over a spherical region from Rmi n , just 
outside the outer boundary of the innermost region with CFR cloudlets, to some 
maximum radius Rmax • The outer boundary, Rmax, may be as large as the radius in 
which the outer thin disk is unstable (see chapter 2). We generally assume ~ 
100/? s , outside of which the thin disk approximation should work (see, for example, 
Frank, King and Raine 1992). Since the IR-optical wave bands have been filtered out 
from the synchro-Compton spectrum by the CFR cloudlets, the radiation pressure on 
the cold clouds due to atomic absorption processes is significantly reduced. Clouds 
larger than the CFR cloudlets can therefore exist in a region larger than Rmin, where 
optically thick clouds are assumed in our work. 

Using the above assumptions the observed spectrum I 0 ba is calculated. The radi- 
ation transfer equation for the modified power-law source around which clouds are 
isotropically distributed (Sivron and Tsuruta 1993). The power-law source can be 
described as /,•„( 0) oc E~ r in most parts of the spectrum, where E = hv is the energy 
of a photon, and T is the photon index. The radiative transfer equation is thus 

Iob,(s) = MO) e~ Tctd + S cld (T' cld )dT' cld , (3.8) 

where s is the path length and 7 in (0) is the input spectrum. The optical depth of this 


assembly of clouds is 


Tdd = f ds'a c id(s'), (3.9) 

Jsq 

which is independent of Tim for opaque clouds. The absorption coefficient for an 
optically thick spherical cloud is a = n c /j7r/^j. The source function Sdd is the energy 
dependent emission coefficient from one cloud divided by its absorption coefficient. 

The outcoming spectrum has three parts. The first part is due to the reflection 
of partially absorbed power law emission from the clouds 

(3.10) 

fir-aicsin(Ftsm9/Rmax) 

exp— / d6 sin ORnddic l m 

JS = o 

where 9 is the angle subtended by the central source and the observer at the cloud. 
The source function S re f(n th (R), I int (R), v) (see, for example, Lightman and White 
1988) is a result of reflection from an individual cloud at R. For example, the solution 
of the radiative transfer equation for low energies is given by S re j = /in(#)(l — 
v /e)/( 1 + y/e), where t = k v I{k u + n ea ) (see appendix 3). The free-free opacity k v 
depends on the density nj tjnH of a cloud, which is derived from equations (2.13) and 
(2.16). of this chapter. 




The additional effects of the relative orientation of the cloud with respect to the 
central point source and the observer are taken into account by the method of moon 
phases (see appendix 2). For example, the decrease in reflection from a single spherical 
cloud is proportional to 7r^jd(l + cos0), similar to the solution for the luminosity of 
the moon phases. The decrease due to this moon phases effect is crudely of the same 
order of magnitude as that of a circular slab tilted at an angle of 9 = 0 — 60 degrees 
(see, for example, the results of the Monte-Carlo simulation in George and Fabian 
1991, Fig 12). The occultation of radiation after it is reflected from the cloud is taken 
into account in the second exponential in equation (3.10). 

The thermal component of the outcoming radiation due to absorption and re- 
emission by the clouds is calculated from: 


1, thermal — f dR [ d0 f 27 rR Sin^ 

JRmin Je = 0 ' *• J 


(3.11) 


{ eip( - L 


?r— arcsin(Hsm B/R m as) 


d6 sin 0Rn c id^l 2 c ii) 


})• 


where F t hcTmai(T c id{R )) is the thermal radiation flux from the surface of the cold 
clouds, and T c id{R) is found from equations (2.11) and (2.19). 


It is convenient to define a covering factor for the central engine. We define the 



covering factor to be 

f(R) = f n c idvlli d dR. (3.12) 

J Rmin 

The covering factor at Rmax is the total optical depth tm for the radially propagating 
photons. If f(Rmax) > 1 our integrals should be expanded to include two reflections 
and more, using a perturbative method. If f{R max ) » 1 the above treatment does 
not work, and the region may be treated as an optically thick photosphere. The 
covering factor at Rmax is a free parameter in our code. 

The observed reprocessed radiation is therefore 

Jobs = 4n(0)e~ T<: ' d + I Te f + I thermal' (3.13) 

The integral of this spectrum, L, is equal to the integral over the primary spectrum, 
because we assumed that there is no internal source of energy for the clouds. The 
emission may be anisotropic if the number of clouds is very small. 

In figure 3-3 we show the radiative transfer in a simplified geometry. The power- 
law source is in the center of a spherical distribution of spherical clouds. The clouds in 
our model absorb, emit and reprocess the power-law radiation. The effects of distance 
from the central source are expected to be implemented through the r dependence of 
the physical parameters calculated in chapter 2. 
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3.2 RESULTS 

In this section we present the spectra from our numerical solution of equation 
(3.13), and explain the consequences of these results. In figure 3-4 we show the 
different spectra for different covering factors. Small covering factor, f(Rm a x ) < 1 
and large optical depth per cloud T\ c \d > 1 are assumed. Note that for f(Rm a x) = 0.9 
the central source is heavily attenuated and the spectrum exhibits many spectral 
features. The spectra in figure 3-4 and figure 3-5 have no R dependence for the 
covering factor, no R dependence for the temperature and R~ 2 dependence for the 
density. The optical depth of each individual cloud is assumed to be much larger 
than unity. We study the conditions under which these assumptions are valid in the 
following two paragraphs. 

The four unknowns n*/,, T c id, n c td and l c u, which determine S re f and F t hermah are 
obtained by solving equations (2.11), ( 2.13), (2.15) and (2.17) for the spherically 
symmetric region. Unlike the case in section 2 of chapter 2 in which we assumed 
/(A*«) << 1 the exponents of equations (2.19) - (2.22) can not be found explicitly. 
It is therefore important to check the limits on the exponents. 

For optically thick clouds the final spectrum is not very sensitive to variations in 
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the radial dependence of the cold density because the 72 dependence of S re j 


is proportional to (1 — \JaRp^~[b + a72 _w )/( 1 + ^JaRp^Jb + aR~ 2S ), where a and b 
are comparable (see appendix 1). The dependence of spectrum on R is extremely 
sensitive to temperature variations because T oc 72\ and the emission is proportional 
to R~ 4X . This result is expected to be limited to the cases of high accretion rates 
in which optical thickness is achieved for most clouds, but may be far from true 
for optically thin clouds with cold matter. The radial dependence of and l c id 
is harder to calculate, but can be estimated with equation (2.17). An upper and 
lower limits on these parameters can be established from as follows: For the upper 
limit we assumed n c id oc 72° and l c u oc 72°. This distribution of clouds, which was 
chosen due to its contrast with the disk model in which no occultation occurs, results 
in a linearly increasing covering factor (equation (3.12)). This 72 dependence might 
be too large far from 72 = R m i n , because the cloud number density must decreases 
with 72 (see equation (2.19)), whereas the cloud size might increase, but then it’s 
opacity tends to vanish as it becomes transparent to the higher energy bands. The 
lower limit is derived by assuming that the equality holds in equation (2.17), a valid 
assumption when the free-fall velocity is proportional to the accretion velocity. The 
radial dependence in the first case is then x = 3, t = 1, £ = 2 and A = 0.5 — 1 and in 
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the second case x = 4, t = 3/2, ^ = 2.5 and A = 0.5 — 1. The second case results in 
a logarithmic increase of f(R) in equation (3.12). The parameter A depends on the 
local covering factor f(R) and is thus not directly related to X an d L - 

The spectra in figure 3-4 thus fall into an intermediate category. (For spectra with 
different R dependence see Sivron and Tsuruta 1994) The spectra were calculated by 
the numerical integrations in equations (3.10) and (3.11). We checked the code with 
four, ten and one thousand iterations. Less then ten iterations are needed in order 
to get the same outcoming spectrum. Because the code is randomly placing clouds 
in 3D space this means that very few clouds are needed per given covering factor for 
there to be isotropic reprocessing, and it would be hard to find the difference between 
few and many clouds. The number of clouds is expect to be larger by a factor of 
l/r XcW for clouds with optical depth smaller than 1. The cloud velocity v cld which is 
uncertain, is assumed to be Keplerian for equation (2.17), an assumption which we 
will justify in chapter §4, but doppler broadening is not calculated (see Sivron and 
Tsuruta 1994). 
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Figure 3.4: The spectra of radiation coming out of a spherical distribution of cool clouds 
are shown with the solid and dot - dashed curves, for the covering factor / = 0.5 and 
/ = 0.9 respectively (the one above). The dotted curve is the primary synchro-Compton 
spectrum. The chain-dashed curve is the spectrum from the innermost region of the central 
engine where very small CFR cloudlets are present (see the text). The dashed curve is the 
spectrum coming out of a slab-corona system. 
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Figure 3.5: The same as figure 3-4, in the intermediate x-ray range. The upper dashed 
curve is for a covering factor of / = 0.01, and the curves below it are for / = 0.5 and 
/ = 0.9 respectively. One can clearly see the ’flattening’ of the power-law spectrum, and the 
increase in equivalent width with increasing /. Also noticeable is the fact that the edge is 
disappearing when / is small. 


The spectrum in figure 3-4 can fit recent observation of radio quiet AGN. These 
observations show a power law spectrum with a bump in the EUV - soft X-rays wave- 
bands, an Fe line edge system and a hard X-ray bump, all of which vary from source 
to source (Mushotzky, Done and Fabian 1993). These results are easily explained by 
a varying covering factor in the cloud model. 

These results, combined with the unabsorbed portion of the input spectrum, can 
be compared with the slab case, which a simplified version of the disk corona model, 
where I obt = 0.5 /,„ + S ref + tt {R 2 max - R 2 min )F th „ mal . In figure 3-4 the dashed line is 
for the slab case. In order to get the same varietyl of spectra the slab must be tilted, 
as in George and Fabian 1990. A disk has a few spectral features which are inherently 
different from those of clouds. The UV thermal bump must be broad, with emission 
that depends on R~ 3 ^ 4 (Frank, King and Raine, 1992), whereas in the cloud model 
the dependence varies from R~ 1 ^ 2 -R~ l . The hard X-ray bump is a result of Compton 
reflection only, whereas in a more general case of the cloud model some clouds have 
r < 1 and a partial covering by a warm absorber should be present (see chapter 7). 
As seen in figure 3-5, the Fe line in the disk model is very strong (except for very high 
inclination angles of the disk, as in George and Fabian 1990). In the cloud model it is 
easy to show that if the Fe line originates in the central engine it may be obscured for 



high accretion rates. This result can be compared with the X-ray Baldween relations 
(Iwasawa et al. 1993), in which higher luminosities were found to be correlated with 


lower line intensities. 


Chapter 4 


CENTRAL ENGINE 
DYNAMICS 


In the first section of this chapter the typical timescales for variations of the phys- 
ical parameters in the central engine are estimated. We show that those timescales 
result in perturbations in the accreting matter distribution, and the formation of 
cold and hot components. In the second section we show that perturbations readily 
form shocks in the central engines of AGN. These shocks are able to reproduce the 
observed flux and spectrum variability, transfer angular momentum away from the 
central engine and accelerate particles to high energies. 

4.1 TYPICAL TIME SCALES 

The dominant processes in the central engine are best determined by the typical 
timescales. Collisional or radiative process are present if their typical timescales 
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are shorter than the typical accretion timescale. The processes with the shortest 
timescales are usually the dominant processes. In this section we try to establish 
timescales and determine the dominant processes. 

4.1.1 Variable Compact Sources 

It is impossible to have highly variable homogeneous compact source. Compact 
sources have very high radiation density which easily accelerate particles to relativis- 
tic speeds. When the radiation is dramatically increased, equilibrium densities and 
pressures are easily perturbed. 

The timescale for doubling the output of the central engine is similar to the ac- 
cretion timescale. If the accretion is homogeneous and the amplitude of variability 
is large enough, the whole emitting region needs to change its emission in a period 
of time t em i 3S i on ~ RJc which is only slightly smaller than the minimum accretion 
time t acc ~ (i2/c)^2r/9, where r = R/R g , without changing the matter density 
distribution. 

The timescale for perturbing the original matter density distribution in a compact 
region is similar to the accretion timescale. One may assume an initial Eddington 
luminosity, such that the matter is at equilibrium. In such a case the doubling of the 
luminosity results in the outwards flow of matter on timescales which are similar to 
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the accretion timescales. Since this timescale is distance dependent the matter closer 
to the source of radiation is accelerated more than the matter further away, and a 
shell, or a sheet, of overly dense matter is formed over timescales which are similar to 
the light crossing timescale. The original homogeneity is therefore easily wiped out. 

4.1.2 Formation and Durability of the Cold Matter Clouds 

In section 2 of chapter 2 the maximum timescale for the formation of cold clouds 
from perturbations in the central engines of AGN is estimated. These clouds should 
persist for a long enough duration so that the effect on the emission is noticeable. We 
next try to estimate the range of lifetimes of the cold clouds using two methods. The 
formation and persistence of cold clouds is then shown to be plausible. 

The minimum life time of plasma sheets and filaments can be approximated for 
the case in which a cloud is, for simplicity, represented by a spherical cloud and all 
of the incident radiation is transformed into an energy increase of size A E. Then, 
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AE 


= 140r'/’ (AltA , 

\ C / XflTUHC*/ 


(4.1) 


where F is the maximum flux incident on the cloud, r» is the radius and T, ~ 10 8 K is 
the final temperature of the cloud in which all the elements are completely ionized (see 
also Bond &; Matsuoka, 1994). For a cloud of size 10 8 cm (obtained for M ~ 10A/6, 
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using the same arguments as in Celotti, Fabian and Rees 1992, or chapter 2 of this 
thesis), the life time of the cloud is greater than 10 -4 seconds. This result is a lower 
limit, because no radiative losses are considered, in spite of the fact that bound- 
free and free-free emissions are dominant mechanisms in the creation of the cold 
component clouds. These clouds have already been proven to be efficiently radiating 
energy in chapter 2. 

The above timescale which describes the final dissipation of a cold cloud is a lower 
limit. We next try to determine an upper limit on the lifetime of a cloud from the 
beginning of such a process. We assume that that the temperature decrease is due to 
adiabatic expansion into vacuum. Since we only have a lower bound on the temper- 
ature of the clouds we get an upper limit on the speed of dissipation. We calculate 
the time-scale for the optical depth to be halved. The following equations govern the 
adiabatic expansion of clouds. The continuity equation in spherical coordinates, for 
a spherically symmetric time-dependent case is 


%^-kh rlpv)=t> ' 


(4.2) 


where r, is the radius of the exploding plasma cloud and v is the velocity. The Euler 


equation is 


dv 1 dP „ 

“77 + _ ”7 — 0 

dt p dr , 

The adiabatic (energy) equation is 

p = Qn r /(i ,-1 ) j 


(4.3) 


(4.4) 


where 7 = 5/3 is roughly the adiabatic index for ideal monoatomic gas, for the mainly 
Hydrogen plasma. The equation of state is 


T ~ 


m H P 
pk B ' 


(4.5) 


The total mass of a cloud, M c ij, is 


A/cW — 



(4.6) 


The optical depth of a cloud, assumed here to be dominated by free-free and bound- 
free absorption processes is 

r ~ 1.0 x ltfVr-'/'zV x (^7) • < 4 - 7 ) 

Through expansion the clouds become more transparent, and their signature on 
the overall emission diminishes. We are therefore interested in the ratio of the optical 
depths of the clouds before and after they expand. By substituting equation (4-8) 
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into (4-7), and than into (4-6) we get T <x r 3 37 , which is identical with the result by 
Zeldovich and Raizer 1967. substituting this result into (4.9) we find: 


iMj o)) _ ( r »(*/) 

T ( r *(*/)) W*o) 


Here rj is the final size of the expanding cloud. From this we can see that a 
bremsstrahlung and bound-free dominated, adiabatically expanding cloud will only 
need to grow by 77 /r 0 ~ 1.12., in order for the optical depth to decrease by a factor 
1/2. Substituting an the initial cloud size into equation (4-5) we can estimate the 
timescale in which a cloud becomes transparent. Since the pressure, density and size 
change very little when a cloud opacity is halved the Euler equation becomes roughly 
cPr^/dt 2 ~ / Po)(P(to) /r*{to)) and integrating twice, using the initial conditions 

v(t = 0) = 0 and r,(i = 0) = r,(0), we get 

t, ~ (lOsec) x (4.9) 

-M> 

where r 8 = r,(0)/10 8 cm, T 6 = r/10 6 K, and r 0 and P 0 are the initial size and pressure 
of the dense sheets, and t 3 is the time from the beginning of the expansion. Hence, 
clouds lose their opacity on a time scale of tens of seconds. The life-time can be longer, 
as a result of confining pressure, due to ambient gas and magnetic fields pressures in 


it (Rees 1987). 
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Optically thick, cold plasmas, form if their creation time is longer than the time 
it takes to radiate the excess energy from the plasmas. In chapter 2 the life time was 
taken to be the infall time. In a more general case, 


< t, 


contraction 


(4.10) 


Here £ e m»jjion> the typical time scale for emission, includes cooling by creation of pair 
cascades or thermal pairs, cooling of the post pair-dominated matter by synchro- 
Compton, bremsstrahlung and free-bound transitions, and thermal emission from 
matter in local thermal equilibrium. The upper limit for the emission timescale is 
given by equation (4-1). Assuming the expansion time is of the same order of magni- 
tude as the contraction time, one can calculate the right hand side of equation (4-12) 
using the expansion timescales t 3 or t m . A minimum cloud size of l c id ~ 0.01T 4 _1 cm 
is needed if t = t ,, and l c id ~ 10 9 r -1 ^ 2 T' 7 -4 cm if t — t *. For t = t s TV = TJ 10 7 is 
the maximal temperature for which equation (4.11) may be applicable. This result 
reaffirms the results obtained by Guilbert and Rees (1988). 

4.1.3 Clouds Collisions 

Clouds do not collapse to form a thin disk only if they are dense and small. The 
timescale for collisions therefore determines the geometrical distribution of clouds. 
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Numerous small and dense clouds may be spherically distributed, whereas large diffuse 
clouds may coalesce and form disks. 

The first simplified approach for calculating the timescale for collisions is pre- 
sented in chapter 2, and is elaborated upon here. In figure 3-3 an ensemble of clouds 
and sheets was presented, to which random motions can be added with an overall 
average angular momentum, with ujR « v avg , where u is the angular velocity and 
v avg is the random speed. The velocity field of inhomogeneities is affected by the 
gravitational and radiative fields, drag by the hot phase, and interactions with other 
inhomogeneities. The average speeds in the case of few inhomogeneities rapidly ap- 
proach Keplerian speeds. Collision between two cold clouds results in the loss of some 
random motion, whereas the overall angular momentum is conserved, and the orbit 
thus becomes more axisymmetric. If the inhomogeneities are not compact the matter 
settles onto an accretion disk. A very stringent upper limit on the number density 
and size of the clouds is therefore found from this argument to be 

(— )”' (*«)■'. (4.11) 

\Vacc' 

Here, unlike equation (2-17), the accretion velocity is not assumed to be proportional 
to vjj. One may assume, for example, an initial distribution of clouds which are held 
at equilibrium (held by radiative pressure, for example), except for one cloud that is 
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only affected by gravity and elastic collisions with identical clouds on the same radial 
line. The result of this may be a ‘domino effect’, in which each cloud transfers all of its 
linear momentum to the cloud next to it, and one gets an equation of motion for the 
radial component 2v avg / At dom = GM/R 2 , where v avg is the average velocity between 
collisions and A t dom is the time between cloud collision in the radial direction. This 
timescale is smaller than A t^u = l/(n c ^7r/^u a „ 3 ), the time between cloud collisions 
in random motion. Comparing these timescales one gets n c i d lh d < l/(2irR), which is 
a stronger limit on the product of cloud density and size than equation (2-17). 

4.1.4 Explanations of Low Amplitude Spectral Variability 

We next show that the lack of simultaneous variability in the IR-optical and the 
X-ray parts of the spectrum does not imply that the synchro-Compton model for 
the central engine fails. In order to show this we establish a relation between the 
temperature of central engine opaque component and the timescale of variability. 

Consider a source in which the amplitude variability is very large in one wave band, 
which includes a significant portion of the total emission. When in certain shorter 
wavelengths the flux variability of time scale A ti ong has small amplitude A Lj At one 
of two things may be the cause: 1. The source of luminosity in that wave band is 
from a region far from the central engine, such that the light crossing time (equation 
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2.1 in chapter 2) is much larger than A t long . 2. The emitting matter has a very high 
‘heat capacity’, such that the energy input from the highly variable central engine is 
reemitted over a large timespan. The mechanism for such ‘smearing’ of the variability 
is described below for the case of opaque clouds. 

For the cloud model a ‘turn-off’ of the central X-ray emitting source results in 
cloud cooling over thermal time scales. The thermal time scale is derived from 
dE/dt = -47r/^ d cr B 7^ (i where E is the total energy of one cloud and cr B = 5.65 x 
10 -5 erg cm 2 K -4 , and T c u(E) is found from the ideal gas law. cr B is the Boltzman 
constant. Solving the above equation we get 

(4.12) 

for Seyfert nuclei, e.g. NGC 4051 (Sivron and Tsuruta 1993). For instance, for very 
soft X-rays at 10 6 K t c ~ 0.01 sec, while it is nearly a year for the UV radiation at 
1000 K. Consequently with our model, while the high energy end of the bump should 
vary almost simultaneously with the harder X-rays, similar rapid variations are not 
expected at the low energy end of the bump lying in the UV band. 

We examined carefully the temporal behavior of the NGC 4051 spectrum near the 
optical B-band (Done et a 1. 1991) which shows no rapid variability during the X-ray 
flares, and find that the low energy end of the UV bump already appears at this 
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tc > 7^3 sec, 



waveband. The lack of rapid variability at these wavelengths is consistent with the 
two component cloud model with no CFR cloudlets if f(Rmax) > 1- This observation 
is consistent with the cloud model with f(Rmax) < 1 if there are CFR cloudlets in 
the central engine, or if the continuum emission mechanism has no inherent IR-UV 
emission (for example, Comptonization of UV photons with relativistic electrons or 
pairs). 

4.2 SHOCKS IN THE CENTRAL ENGINE 

Shocks form readily in the central engines of AGN. Shocks are shown to form 
when a star is captured into the central engine, or directly, from strong and large 
density perturbations. If the shocks persist for a long enough period they are shown 
to significantly modify the emission from the central engine. Shocks also provide the 
necessary mechanism for the outward transfer of angular momentum which enables 
the accretion of matter onto the SMBH. 

4.2.1 Introduction to Shocks 

Using recent X-ray observations it was shown that a large fraction of the high en- 
ergy luminosity is emitted by hot ( T e > 10 9 K) plasmas (Blandford 1990, Mushotzky, 
Done and Pounds 1993, Svennson and Zdziarsky 1994). In order to explain the in- 
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creased emission from these plasmas a triggering mechanism should be devised in 
which a small perturbation can generate radiation from a large enough volume. The 
effect of the perturbation on the hot plasmas within this volume double its’ radiative 
output, in some cases. 

The largest speed in which perturbations can propagate is of the order of the speed 
of sound or the Alfven speed in the hot plasmas. Perturbations that are large enough 
to result in significant changes in emission are easier to produce by shocks than with 
regular compression waves. We therefore look for a mechanism that results in shocks. 
One should check that shocks indeed form under the conditions anticipated in the 
central engine, which we do in section 4.2.2. We show that shocks form in the central 
engines of AGN. We discuss what happens when a local cold inhomogeneity, such as 
a cloud, or a main sequence star, makes its way into the innermost region of AGN. 
We have to show that such an incident naturally results in the creation of a strong 
shock that may modify the emission from the continuum source, the hot plasma. We 
have chosen a main sequence star for two reasons. The first is that the capturing of 
such a star is a plausible event (Syer Clark and Rees 1990), and the second is that 
an analytic solution is possible when the lifetime of the inhomogeneity is long, as in 
the case of a star. We also discuss the consequences that the shock might have on 
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the accreting matter and the star. A more direct approach, in which a non-linear 
perturbation results in a shock wave, is calculated in our next paper (Sivron and 
Tsuruta 1995). 

Finally we have to show that a relatively small inhomogeneity significantly mod- 
ifies the light-curve from an X-ray region of at least an order of magnitude larger 
surface area. The creation, destruction and properties of the emission modifying 
matter are treated in section 4.2.3 The resultant flux variability is shown to be sig- 
nificant in section 4.2.4 A simplified geometrical model used to calculate a typical 
light-curve is presented in section 4.2.5. The problem of angular momentum transfer 
is treated in section 4.2.6. Shocks, turbulence and magnetic reconnections are prob- 
ably the primary mechanisms in AGN for the creation of cold clouds, as is further 
discussed in chapter 7. 

4.2.2 Shock Formation and Its Consequences 

We next show that the velocity of a macroscopic inhomogeneity exceed the speed 
of sound and/or the Alfven speed when captured into the hot plasma in the inner 
region of an AGN. A shock wave, which always forms when matter supersonically 
flows past a body (Landau and Lifshitz 1987, hereafter referred to as LL87), forms 
near the inhomogeneity. A collisionless shock will form if the Alfven speed is exceeded. 
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Theory predicts (Shapiro Lightman and Eardley 1976, Rees et al. 1982, Blandford 
1990) that the thin accretion disk will be disrupted within a critical radius, which may 
be as small as r CT ~ or as large as r^ = m p /me, where a is the 

parameter defined in the ‘standard a disk’ model, which relates the pressure and shear 
stress. Within r cr heating due to gravitational and viscous energy release surpasses 
cooling by bremsstrahlung and black-body radiation, the temperature rises and the 
gas pressure inflates the disk. The accretion flow inside this radius will therefore 
form a geometrically thick torus, which will become a two temperature optically 
thin plasma for rates below the critical rate rha- ~ 50 a 2 (Rees et al. 1982). Even 
at higher accretion rates a two temperature state may ensue (White and Lightman 
1990), which is not necessarily stable (Begelman, Sikora and Rees 1987, White and 
Lightman 1990). We hereafter concentrate on the optically thin two temperature 
plasma, and use the equations of Shapiro Lightman and Eardley 1976 and White 
and Lightman 1990. The thermodynamic parameters, such as the adiabatic constant 
T a , are taken from a paper on shocks in pair-rich plasmas (Iwamoto 1989, hereafter 
referred to as 189 and references therein). Although the 189 equations, which are for 
matter in LTE, may be over-simplified for parts of the central engines of AGN, the 
results are expected to be qualitatively correct. 
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The speed of sound at a distance r < r c- from the black hole is: 


c, = (4.13) 

V 

where p is the pressure in the accreting gas, and u; is the enthalpy he 189 equation 
A18a). The pressure is dominated by the ion thermal energy, and the equilibrium 
temperature does not quite reach its virial value if the outwards transfer of angular 
momentum is dissipative. (This is because of a combination of inverse Compton 
scattering and bremsstrahlung in mild pair plasmas, synchro- Compton cooling in 
strong magnetic fields, or some collective plasma effects, such as those created by 
shocks.) The ion temperature is thus 

Ti = a T v , (4.14) 

with a ~< 1 and T v = m v c 2 /(fear) is the virial temperature. The ideal gas law, which 
holds for optically thin gasses at high temperatures (189 eq. A8), is a reasonable 
approximation for the equation of state. The speed of sound is therefore 
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where mjj is the mass of the hydrogen atom and for cosmic abundance /x ~ 1. 

An inhomogeneity, which is more compact than the surrounding material, so that 
it cannot lose its angular momentum on dynamical time-scales (see sections 4.1.2 and 
4.1.3), moves at a speed comparable to its Keplerian or free fall velocity, both of 
which are usually larger than the speed of sound: 


t>, - |v. - v,/| = 4r 1/3 c. 


(4.16) 


Here v* is the velocity of the inhomogeneity and ~v fi is the velocity of the accreting 


matter. The Mach number 




(4.17) 


does not depend on the distance from the SMBH, because both speed of sound and 
typical velocity are proportional to the virial energy. For b ~ 5/4 the Mach number 
is larger then 1, and a shock forms. Because a is not much smaller than 1 for a 
wide range of accretion rates (see WL90) M is not larger than ~ 10, and a shock, 
rather than tube, form. The Mach surface opening angle, <j> ~ 2 arcsin(c,/u*:), is 
therefore significant. The shocked surface so created roughly takes a conical shape. 
Higher accretion rates will make Y closer to its radiation dominated value of 4/3, 
which makes shocks possible even in the cases of near virial temperatures and average 
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inhomogeneity velocity. 

The creation of strong shocks by large enough inhomogeneities in the central 
engines of AGN is inevitable for a large range of accretion and dynamical parameters. 
These inhomogeneities can be stars, but may also be cold plasma sheets and filaments, 
that accelerate to supersonic velocities if the gravitational force, the radiation pressure 
or the (virial) magnetic stresses are greater then the drag force 


GMm* 

B? 


> irBjpo, 


(4.18) 


where p 0 is the pressure on the inhomogeneity surface (see section 4.2.3). We assumed 
that the radiative force or the magnetic stresses are of the same order of magnitude 
as the gravitational force (see also chapter 2). Noting that the pressure is roughly 
po ~ Pg v K ( see section 4.2.3) this condition can be also written in the form: 


->-£ (— ), 

Pg D.\p) 


(4.19) 


where p, is the average density of the inhomogeneity. This condition, which is fulfilled 
for inhomogeneities which are confined by virial magnetic fields (see Guilbert and 
Rees 1988, Sivron and Tsuruta 1993, and Celotti, Fabian and Rees 1992), may be 
less restrictive if we include acceleration by radiation. If confined by virial magnetic 
fields the inhomogeneities may be as small as D+/Rx = 10~ 4 , two orders of magnitude 
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smaller than an orbiting star. The amplitude of perturbations in the density needed 
for shocks to form is therefore small, and shocks may form in the wake of large 
amplitude compression waves in the disk. These waves are probably non-axisymmetric 
(see section 4.2.6). 

It is easier to calculate the resultant shocks for inhomogeneities with lifetime 
which is much longer than the dynamical timescales, such as stars. It has been shown 
that a star can be captured into a tightly bound orbit around a black hole through 
repeated interactions with a thin disk around the hole (Syer Clarke and Rees 1990), 
and that the probability of such an event in AGN is even higher for stars in bound 
orbits (Pineault and Landry 1994). Once captured the orbit will be circularized and 
eventually the star will be “grounded down” into the disk plane through repeated 
interactions with the disk. The timescale for circularization can be shorter than that 
for planarization, which may happen very close to the central source. Therefore, 
in our model, we assume that the orbit of the captured star has been circularized, 
whereas the orbital plane is still inclined to the disk plane. This situation is expected 
for a wide range of the original relative inclination angles of these planes. 

The extent and magnitude of the shock depend on the size of the star, and so does 
the stability of the orbit. For example, for a similar orbital velocity a smaller star 
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produces weaker shocks than a larger star, the drag force is smaller, and the decay of 
the orbit is slower. To avoid rapid disruption the size of the star that has migrated 
into the innermost region cannot be much larger than its Roche lobe. The mass of a 
regular main sequence star is therefore limited to 

m* < 0.3 Mq (4.20) 


(King & Done 1993). The radius of the star may therefore be expressed in units of 
d\o = D»/10 xo , and its mass is in units of the solar mass Mq. 

In most cases the generalized hydromagnetic speed of sound, ( c*0 ) 2 = cj + v\, 
where va = 8xp) is the Alfven velocity and B is the magnetic field, is similar 

to the hydrodynamic speed of sound (equation (.•••))• The main effect that magnetic 
fields have is reducing the shock thickness: 



(4.21) 


where a ~ ne 2 tc/m e , tc is the ion-electron collision timescale and B\-i is the mag- 
netic field variance across the shock. The smallest possible shock width is the Larmour 
radius, which for ions is ru = AT x ^ 2 jZB ~ \Q^Tll 2 j Bern where Am p is the mass of 
an ion and Z is its charge. It is evident that even magnetic fields which axe much 
smaller than the virial value (about 10 3 Gauss in the central engine, see Bland ford 
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1990) limit the thickness of shocks. 

Although the speed of sound in The extent and magnitude of the regular super- 
sonic shocks depend on the size of the perturbation and the density of the perturbed 
matter, and in the case of AGN intermediate and high accretion rates make creation 
of significant shocks possible only in the inner 1072 g- In that region the mean free 
path in the surrounding gas is smaller than the size of the inhomogeneities. How- 
ever, if the accretion is not super Eddington, collisionless shocks form throughout the 
central engine, because the Alfven speed is usually much smaller than the speed of 
sound. 

The Alfven wave speed in plasmas is roughly C\ = yjB*/8Tcn 3 m a , where the 
subscript s is for electrons or ions. Super Alfven speeds are therefore achieved if 
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is larger than 1. Here nio = n 4 /10 lo and B3 = B/10 3 . For the central engine one 
expects B fields of order 10 3_4 gauss or smaller (Blandford 1990). It is apparent 
that the Alfven speeds are exceeded in most cases. Large magnetic fields and small 
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densities may be present in radio loud AGN. In such cases the ordered structure 
of the accretion disk is not perturbed by shocks, and advection, magnetized disks 
and jets may be possible. For higher accretion rates diffusion of magnetic fields out 
of the plasma and reconnection events may restrict the magnitude of these fields. 
We therefore expect that for near-Eddington accretion rates the magnetic fields are 
smaller than 10 3 gauss, and collisionless and magnetic shocks form. 

4.2.3 Radiation From the Shocked Matter 

We show that the when the power transferred from the star to the shocked matter 
is emitted as radiation the light curve is modified. We then show that its highly 
probable that the shocked matter indeed radiates the power transferred from the 
star. 

The power transferred from the star to the shock is at least 

Pst s ~ (A,po) x u„, (4-24) 

where A, = nD\ is the cross section area of the star and po is the pressure near the 
stellar surface, behind the shock. The pressure on the star surface, po, is modified 
from the treatment of formation of shocks in supersonic non-relativistic flows past 
solids of revolution (LL87, §122). In our case the object is a sphere, the orbiting star. 



The pressure near the stellar surface is thus related to the ambient matter velocity, 


density and speed of sound through 

»p a _j_ j\ (r<.+i)/(r«-i) 
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where p g = M/(QR 2 v acc ) for semi - spherical accretion, ft is the solid angle and v acc 
is the accretion velocity. If the parametric dependence of the density is less distance- 
dependent (see, for example, WL90) the pressure is even higher at larger radii. The 
adiabatic constant T was chosen to be at it’s highest value of 5/3, which results 
in a lower limit on p 0 . We have used the Ranking-Hugoniot relations for relativistic 
plasmas in thermal equilibrium (Iwamoto 1989). For collisionless shocks the minimum 
pressure on the stellar surface is po = B 2 /^x. 

The power transferred from the star to the shock is thus at least 

P»u ~ ( 10 41 er^s sec - 1 ) 6 3 d\ Q M 2 3 Me 2 r~ 3 . ( 4 . 26 ) 


If all this power is radiated the radiation emitted is quite large. This radiation may 
be smaller than the X-ray luminosity observed (table 1), but it may still modify the 
light curve and spectrum. We also note that relativistic effects, secondary shocks (see 
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section 4.2.4) and the effects of short living larger inhomogeneities are not taken into 
account. The modification of the emission through a single radiative shock may occur 
only in the innermost 10/? s , as is anticipated (but not calculated) by previous authors 
that discussed pair modifications of the inner torus with the injection of a power-law 
distribution of relativistic electrons. (See Blandford 1990, White and Lightman 1990, 
Tritz and Tsuruta 1989, Kusunose and Takahara 1989, Kusunose et al. 1995. These 
authors mistakingly assumed strong shocks exist in the central engine by a-priory 
postulating geometrically thin, optically thick Keplerian accretion disk. In such disks 
the speed of sound is much smaller than that calculated in equation (4.12), but these 
disks are not necessarily stable.) 

To show that the power transferred to the shock is indeed radiated one needs to 
look at the Ranking-Hugoniot relations away from the star. The increased densities, 
temperatures and particle number densities in the post-shock region must be large 
enough so that the timescale for radiative processes is smaller than the timescale for 
the dissipation of these quantities. For strong MHD shocks the mean free path for a 
particle in the post-shock region is approximately the shock width A au , ~ 10 11 nf 3 1 cm 
(see standard plasma texts such as Boyd and Sanderson 1968), which means that the 
radiative cooling timescale must be shorter than the X 3W /v m ~ lOOnj^r 1 / 2 seconds in 
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which the matter passes in the post-shock region. In Guilbert and Rees (1988) and 
Ferland and Rees (1988) the cooling of a slightly compressed region in the central 
engine down to LTE temperatures by bremsstrahlung and synchro-Compton radiation 
was calculated to be possible for a region which is two orders of magnitude larger. 
Their calculation therefore assumed the average densities which were calculated in 
chapter 1. 

In order to calculate the Ranking-Hugoniot relations in the shock we assumed 
thermal two temperature plasmas. This is a reasonable assumption for high accre- 
tion rates and moderate magnetic fields. This kind of formalism was developed by 
Iwamoto (1989) for plasmas near pair-creation temperatures, which is appropriate for 
the central engines of AGN. Using this formalism for the post-shock region the density 
is calculated from the Ranking-Hugoniot relations to be of order 7 p g . The post-shock 
ion temperatures, including energy deposited into pairs created in the shock, are of 
order 107) for the central engine conditions of table 1. The synchro- Compton cooling 
(and pair cooling for high compactness sources) time is roughly shortened by at least 
taw/tGR rsj where the subscripts GR stand for the Guilbert and Rees 

result. The timescale t aw is the time for cooling in the shock. The factor 1/7 comes 
from the increased post-shock density effects on the synchrotron and Comptonized 
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radiation. The temperature of a shocked particle T oc KE ~ mc J (7i — 1 ) implies 
that the Lorenz factor for the average downstream particle, 7 pj , is roughly an order 
of magnitude larger than the upstream 7 gr- It is thus evident that strong shocks 
can result in cooling timescales that are shorter than the typical crossing time of the 
shocks. The matter that was cooled by synchro-Compton processes can be further 
cooled by bremsstrahlung radiation. The timescale of cooling by bremsstrahlung is 
proportional to the square of the density, and is therefore much shorter than the time 
for reheating by radiation (see chapter 2). This cool matter is probably sheared and 
fragmented in the rarefaction region behind the shock. We see that not only does 
the post-shock matter can radiate most of the energy transferred from the star to the 
shock, but also post shock thermal (and possibly opaque) Guilbert and Rees ‘cold’ 
matter may forms. 

In the above treatment we assumed that the magnetic field is perpendicular to 
the shock front and is unchanged in the shock. In such cases the shock thickness S 
equals the mean free path of the ions (or the deflection radius) rather than the much 
smaller Larmour radius in the magnetic field. 

In order for the power from the star to be given to radiation from the hot shocked 
matter the mean free path must be small enough. In the case of collision shocks the 
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mean free path is larger than the star, and smaller than the emitting region. In such 
cases shocks are inefficient. However we argue that magnetic fields can be rid off only 
to the level at which they are reconnected in reconnection events, such as in shocks. 
Since magnetic fields can be refracted and changed by shocks current sheets and DC 
fields can form. These can be the sight of reconnection events. If such is the case one 
expects that the leftover magnetic fields are enough for Larmour radiuses which are 
as wide as the typical inhomogeneity, which we assumed to be of order 10 lo cm. The 
timescale for cooling is still enough, because the magnetic fields are increased by the 
compression ratio rc = p?/ P\ ~ 7. 

The preferred shock thickness is therefore much smaller than the mean free path. 
Indeed it may be as small as the Larmour radius of an ion in the field. The case 
of collisionless shocks, which are expected in AGN for low accretion rates, is slightly 
more complicated because of the supposedly small shock front width and the increased 
importance of Fermi-shock processes, and because of the breakdown of the thermal 
plasma approximation (see Sivron and Tsuruta 1995). 

The shock induced modifications of the emission from the central engine may 
result in a flare of softer than synchrotron X-ray spectrum, i.e. 0.5 < A < 1, where 
A is the slope of the continuum emission spectrum, because of the increased local 
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compactness (see, for example, Ghisellini and Haardt 1994, Grandi, Done and Urry 
1994, and references therein), or because of the spectrum from reconnection events. 
Temporary increase in the brightness of an inhomogeneity could also be a result of 
lensing of the shock radiation near a Kerr black-hole (Blandford and Narayan 1993, 
Rauch and Blandford, 1993). 

Significant amounts of radiation from the shock front will also be observed if more 
efficient non-thermal pair cooling is present. Such cooling will result, for example, 
if the Fermi acceleration is very efficient (Blandford and Eichler, 1987, show that 
electrons heating by interaction with Alfven waves is limited, but relativistic pairs 
from various pair-creation events can efficiently cool the post shock matter by synchro- 
Compton and bremsstrahlung processes). 

4.2.4 Effects of Post-Shock Matter on Radiation 

If P ata is expended on the cooling of the post-shock matter, that cool component 
may modify the emission by absorption in a post-shock sheet behind the Mach surface. 
We next show that the power needed to cool a large enough sheet is smaller than P ata . 
The X-ray continuum source, usually assumed to be ~ 10 R g in radius, 


R x ~ (1.5 x 10 12 cm)M 6 , 


(4.27) 
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should be covered by a large enough sheet of area larger than about nfi? x . The 
continuum region can be obscured if the absorption column through this sheet is 
larger than Njj ~ 10 2O-24 cm -2 , where the lower number is for bound-free absorption 
in cold (the ionization parameter Ho = 0) cosmic abundance matter in which some 
features, such as an oxygen and iron K edges and lines may show up, and the higher 
number is for electron scattering. The sheet thus created must therefore be made out 
of at least NhkR x = 7 x 10 44_48 M| particles. 

In order that pair cascades will efficiently cool enough of the ambient matter the 
energy per particle, e ~ (3/2)£flT 2 should be ~ 10 m K <? (see, for example, Svennson 
1984, Lightman 1982), where T 2 is the post shock matter temperature. The product 
of the rate at which particles are cooled and the life time of the post-shock cooled 
matter should be larger than Nu. (The life time of a cloud, t 3 , is calculated in the 
Appendix.) The power of the shock therefore supplies enough energy per particle for 
cooling to take place 

v Ratals 

= (0.1 — 1000)m e c 2 ^-d 2 0 M 6 ' 4 M 23 r -3 . 


(4.28) 


(The above quantity is similar to the ‘compactness’, l = LaTl(Rm e <?), which was 
defined by Guilbert Fabian and Rees 1983, where cry is the Thompson cross section. 
One may thus think of shocks as an efficient mechanism for the local amplification of 
the compactness.) If we consider the synchrotron power P ayn 1° equation 4.25 instead 
of P aia we will still get large enough energy per particle if the post shock magnetic 
fields are increased by a factor, as is expected in collisionless shocks. 

Cold matter which is optically thick in soft X-rays may therefore exist in regions 
which are inside the thick disk if the inhomogeneity is as large as a star. (It is much 
more likely, though, that this matter will be ‘warm’, i.e. partially ionized by the 
stronger radiation from the central disk and from the nearby shock.) In the case of 
nonthermal pair cascade the existence of this matter is even more plausible. This 
cooled matter can be the same as the ‘cold phase matter’ suggested by Guilbert and 
Rees 1988 (hereafter referred to as GR88), and may take a spiral shape (see discussion 
in section §5). We note that modifications due to Fermi acceleration, will probably 
lower the efficiency of cooling while enhancing particle acceleration for the higher 
energies in equation (11) (Blandford and Eichler 1987, Levinson 1994). 

4.2.5 The Effects of Shocks on the Overall Flux Variability 


Our result for the minimum power transferred by the star to the shock may be 
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used for the purpose of finding an overall flux variability. The simplest scenario is 
thus a slow introduction of ‘stars’ of different sizes (see also section .... for details on 
the source for such ‘stars’) into the central engine with no angular momentum. The 
flux from each such event is P a t,tin which can be written in the form 

U = (3 x 10 42 er^)6 3 d? o M 23 M 6 - 2 r- 3 , (4.29) 

where rj 2 is the distance from the SMBH in units of 10 12 cm. A random intro- 
duction of cold inhomogeneities with random distribution of Areas A = 10 2O 7rd? o 
naturally results in shot noise, providing that the events are not correlated. 

When the accretion rate is higher the events are correlated, for example, by means 
of radiative pressure and radiative transfer. The speed of the infalling cold inhomo- 
geneity is reduced, lowering P ata , while the life time is reduced because of heating 
which speeds the initially adiabatic expansion, thus reducing the length of period of 
the event, and limiting the amplitude of variation. This scenario is similar to self or- 
ganized criticality in phase transition physics (Bak et al. 1989). The power spectrum 
from such an event is expected to be a chaotic in nature, i.e. the slope of the power 


spectrum will reach — 1. 



4.3 OCCULTATION BY POST-SHOCK MAT- 
TER 

In this section we calculate the effect of an optically thick sheet behind a shock 
on the light-curve. The main purpose of this exercise is to demonstrate the type of 
spectral variability to be expected when shocked matter modifies the emission from 
the central continuum source. Our approach is geometrical in nature. In our model 
the body obscuring the source is not circular. We treat the cases in which a triangular 
shock front eclipses a circular or a square continuum source. 

Possible picture of what happens near the central region of an AGN with an 
inhomogeneity is shown in Figure 4-1, in which we see the shock inducing star (or 
spherical inhomogeneity), at point S, orbiting inside the thick accretion disk outside 
the X-ray continuum region, but confined to a region smaller than the thin disk. 
From an observer’s point of view the points A, S and B on the Mach cone form a 
‘fan’ in the shape of a triangle that partially obscures the internal continuum source. 
The type of occultation depends on the optical depth of the matter on and near the 
shocked surface. Emission from the shock surface may dominate the light curve in 


some configurations. 
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Since the extent of the optically thick shock may be limited, as was shown in a 
previous section, we choose the point of view of an observer in direction G (Figure 
la. and lc.), which requires the shortest extent of optically thick matter. Figure lb. 
is applicable only when the shock effects extend far enough from the star. 

Figure 4-2 is a simplification of the view from the direction of point G which is 
used for the calculation of the light curve. It is based on the following assumptions: 

1. The continuum X-ray emitting region is assumed to be circular, with homogeneous 
surface brightness, as viewed when unocculted. Its radius is R — Rx • We note that 
in an alternative model, in which the soft X-ray region is larger than the hard X-ray 
region, the soft minimum flux will be significantly delayed with respect to the hard 
flux. 

2. The Mach surface is assumed to be an equilateral triangle, passing across the 
center of the circular source. That is how a conical surface looks to an observer in a 
direction perpendicular to its altitude. 

3. Density variations are simplified by assuming three density regions: in region 
1 (angles smaller than (f>\) the matter is opaque to both hard and soft X-rays. In 
region 2, with angles between <f>\ and <f > , the matter is transparent to hard X-rays, 
but opaque to soft X-rays. In region 3 the density of matter varies with distance, 
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due to the spread of matter, and is partially transparent. We consider only the hard 
X-ray wave-band in the 10-20 keV window in which no absorption is expected, and 
the soft X-ray wave-band in the 0.5-2 keV window in which absorption due to cold 
or warm oxygen is expected. We chose the maximum covering factor to be / = 0.5 
for demonstrative purposes. 

In Figure 4-3 we see the light curves near a dip in the two energy bands that will 
be observed when the source is eclipsed by the post-shock matter. The softer (dashed 
dotted) and the harder (solid curve) X-ray windows with <j> = 45°. The dotted curve 
is for the case of a smaller Mach surface opening angle of <f> = 30°, for the hard X-rays. 
The dashed curve is for the simplified case in which a triangle with an opening angle 
<j> — 45° eclipses a square source. In all cases the maximum covering factor is / = 0.5. 
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Figure 4.2: An observer’s view of the AGN in our model with covering factor of / = 0.5. 
This is a simplified form of the view from point G in Figure 1. The stepwise fashion which 
was chosen to describe the absorber is for demonstration purposes. See section §4 for details 

c 



Figure 4.3: The light curves near a dip that will be observed when the source is eclipsed by 
the post-shock matter is shown for the softer (dashed dotted) and the harder (solid curve) 
X-ray windows with <f> = 45°. The dotted curve is for the case of a smaller Mach surface 
opening angle of <p = 30°, for the hard X-rays. The dashed curve is for the simplified case 
in which a triangle with an opening angle <f> = 45° eclipses a square source. In all cases the 
maximum covering factor is / = 0.5 (See section §4 for details). 
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The decrease of flux toward the dip is energy dependent. A delay is observed 
in the sense that the decrease in photon count in higher energies lags the decrease 
in the lower energies. The light curve in the hardest window is chosen to have a 
cusp-like shape, dropping to half its value at the bottom so that the size of the 
central sector of absorbing material within <f>i which is opaque to the harder X- 
rays is approximately half the size of the harder X-ray source size during maximal 
occultation. The maximum value for Rx is thus obtained with <f> i = tt/ 4. A similar 
method is used for the simpler geometrical model. It is therefore speculated that the 
spectral variability observed in AGN may be with time scales which are related to 
the distance of post shock matter from the central continuum source. 

4.4 TRANSFER OF ANGULAR MOMENTUM 

In AGN shocks seem to be an efficient mechanism for cooling (but shocks in 
standard thin disks cannot self consistently explain the high-energy spectrum, see 
Chakrabarti and Wiita 1992 and references therein). The amplitude of perturbations 
needed for shocks to form is therefore rather small, and shocks may form in the wake 
of large amplitude non-axisymmetric waves in the disk. 

Shocks may thus be also an efficient mechanism for the transfer of angular momen- 
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turn: Spiral shocks may form when a massive companion perturbs the disk (Morfill, 
Spruit and Levy, 1993, and references therein, in reference to galactic sources). Such 
companion cluster which was predicted to exist near AGN (Rees 1990 and references 
therein), was recently observed in M31 by the Hubble Space Telescope (HST) (Pei et 
al. 1994). Such shocks are considered to be an efficient mechanism for the outwards 
transfer of angular momentum in the case of galactic and proto-stellar disks when the 
mach number is not too large (A4 < 10, Savonije, Papaloizou and Lin, 1994, Morfill, 
Spruit and Levy, 1993). This possibility should probably be reconsidered for AGN 
due to the HST observation of a spiral structure in M87 (Ford et al. 1994). It must be 
noted, however, that spiral shocks will result in non-chaotic X-ray temporal variabil- 
ity in low accretion rate AGN, that may be revealed by long exposures observations 
by X-ray telescopes such as Rosat, Asca, XTE, XMM and especially AXAF. The 
angular momentum transfer in the central engine in high accretion Seyfert 1 galaxies 
and radio quiet quasars can be due to chaotic, rather than spiral, shock structure. 
(In these sources temporal variability is already known to be nearly chaotic.) 

A comprehensive study of the transfer of angular momentum in the outer disks 
(Caditz, Sivron and Tsuruta 1995) and central engines (Sivron, Caditz and Tsuruta 
1995) with shocks is currently under way. 
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Chapter 5 

OBSERVATIONAL METHODS 

In this chapter we describe some observational methods in X-ray astronomy. Since 
we are interested in the analysis of data from X-ray satellites we emphasize the pos- 
sible systematic and random errors in the instruments on board the satellites. The 
observation of weak X-ray sources, such as AGN, also requires the development of 
careful methods for the analysis of the data. We describe those methods for the Ginga 
and ROSAT satellites. 

5.1 X-RAY DETECTORS IN ASTRONOMY 

In this section we describe the physics of the proportional counter, the ‘work horse’ 
of X-ray astronomy. By understanding proportional counters one can understand 
most of the errors and background subtraction in X-ray data. (In this thesis only 
data from proportional counters is analyzed. Other X-ray detectors and techniques 
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are described in Fraser 1988.) A brief description of the detection of X-rays is outlined. 

X-rays in the range 0.1 — 20keV entering a detector are absorbed by electrons 
in the innermost shell of an atom with high enough Z. Atoms with large enough Z 
have binding energies in the X-ray range for their innermost shells. Noble gasses, 
Argon and Xenon, are used so that chemical reactions, which complicate the energy 
transfer, will not be prevalent. The energy of the incoming X-ray photons is usually 
transferred to a cloud of ions and electrons through a combination of photoelectric 
and Auger processes, and later they are multiplied by ionization avalanches in the 
strong electric field around the anode wires. For example, a photo electron is ejected 
from the inner most K-shell. The vacancy in the K-shell thus created is filled by an 
electron in outer shells. The energy difference between the two shells is released in one 
of two ways: Auger electron emission from the outer most shells, or fluorescent X-ray 
emission. The primary photo electrons and Auger electrons dissipate their energy by 
ionizing electrons in outer shells of other counting gas atoms. Finally, all the energy 
of the incident X-rays, E(KeV), is converted to a primary electron cloud with total 
number N, which is proportional to the primary X-ray energy as follows, 


" = # 


(5.1) 


where W is the mean ionization energy of the counting gas, 26.2eV for Argon and 
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21.5eV for Xenon. Along the electric field of the detector, these electrons are collected 
and accelerated toward the anode wire. When the electron obtains enough energy 
to ionize another counter gas atom electron multiplication process takes place (this 
happens within a mean free path distance). Usually one primary electron yields on 
the order of ten thousand electrons on the anode wire. Therefore the height of the 
pulse is proportional to N x M, where M is the gas multiplication factor. Such 
proportionality of the pulse hight to the incident X-ray energy is the primary reason 
for the name ‘proportional counter’. The energy resolution of about 18% results from 
fluctuations in the primary electrons number and the gas multiplication factor. 

The detected energy, number of photons and the time of detection are recorded by 
the satellite and relayed by telemetry to the relay station on earth. The information 
is stored on magnetic tapes, usually in a binary form. 

5.2 THE GINGA MISSION 

In this section we describe the Japanese X-ray satellite Ginga (1987 - 1991). The 
large area counter (LAC) on board Ginga, which used multi wire proportional counters 
technology, is described. Background rejection methods for LAC are explained, and 
techniques of data analysis are outlined. A more detailed technical description of 



Ginga can be found in Turner et al. 1989 and Makino 1987. Hayashida et al. 1989 
and Awaki et al. 1991 explain the background rejection methods for Ginga in detail. 

5.2.1 Satellite Description 

In figure 5-1 we see the satellite Ginga. 


All Sky Monitor 



Figure 5.1: The Ginga satellite: the LAC eight detectors are pointed towards the left. 
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Ginga was launched in February 1987 into a low earth orbit with perigee of 505 
km and apogee of 673.5 km, inclination of 31.1° to earth’s longitudes, and orbital 
period of 96.5 minutes, which is an equivalent of roughly 15 orbits per day. Ginga 
satellite orbit decayed slowly as a result of interaction with the upper atmosphere 
(drag), and was lost on October 1991 when the increased solar activity resulted in 
increased drag. 

The characteristics of Ginga’s orbit determine its observational capacity. The 
relatively low orbit protected the satellite from earth’s radiation belts, except in the 
South Atlantic Anomaly (SAA) where the radiation belt is ‘hanging down’ because 
of the weak geomagnetic field. Because of the low orbit the satellite body was also 
protected from most high energy cosmic rays which are deflected by the geomagnetic 
field, below energy referred to as the Cut Off Rigidity (COR). Regions with small 
COR are therefore subjected to high particle background. Because of the high particle 
background in the SAA weak pointed sources, like Seyfert galaxies and QSOs, could 
not be observed in the 5 orbits that transcended the SAA. This limited the observation 
of AGN to 10 orbits per day. In low earth orbit all sources are observed intermittently, 
because of earth occultation of the source. Out of the 96.5 minute Ginga orbit up to 
~ 65 minutes could be utilized for pointed observation, depending on the coordinates 



of the observed source. During the rest of the period view of the source is blocked. 

the three astronomical instruments on board Ginga were the large area counter 
(LAC), the all sky X-ray monitor (ASM) and the Gamma-ray burst detector (GBD). 
Here we present a very basic description of LAC (for further details see Turner et al 
1989), which is the only Ginga instrument used for pointed observations of AGN. LAC 
consisted of eight proportional counters arranged in two layers with total effective area 
of 4000 cm 2 . The counters were filled with a mixture of 75% argon, 20% xenon and 
5% CO 2 . The x-rays entered through a 62/xm thick beryllium window which absorbed 
very low energy x-rays. The top layer was sensitive in the 1.5 — 37keV range, while 
the ‘MID’ layer was sensitive in the 5 — 37keV range. The combined sensitivity of 
the two layers is shown in figure 5-2. LAC is efficient in the 2 — 20keV for weak 
sources. The angular field of view with the detector was 1.1° x 2° degrees, with the 
long axis perpendicular to the satellite equator. LAC had 64 energy channels, out of 
which the two lowest energy and the upper channels were not used for pointed source 


observations. 
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Figure 5.2: The relative effective area per photon energy for a combination of top and 
layers LAC on Ginga. 
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LAC on Ginga was a very successful instrument. The LAC detector had the largest 

effective area of any detector that has flown so far in this energy range. As such, it 

collected a high number of photons and could determine flux level, spectral features, 

flux variability and spectral variability with high statistical significance. LAC was 

therefore useful in discovering new spectral features in AGN. It was also the first X-ray 

detector that had significant sensitivity in the 10 — 20keV energy range (Mushotzky, 

\ 

Done and Pounds 1993). Data from LAC was also used to determine significant X-ray 
spectral variability in AGN for the first time. 

5.2.2 Background Rejection Methods 

Unlike telescopes a large area counter has similar collection and detector areas. 
This makes background rejection very complicated. The many events detected by the 
large area detectors which are the result of ionization of and decay of detector parts 
by high energy particles and direct detection of particles is referred to as internal 
background. The large angular field of view requires a very accurate subtraction of 
the X-ray background radiation and radiation from other sources in the field of view, 
referred to as external background, which is not spatially smooth on this angular 
scale. Ginga background subtraction is especially challenging, because no simultane- 
ous background observation was performed. For weak sources the number of total 


events exceeds the number of source events even after sophisticated anti-coincidence 
techniques are applied. A complete description of background rejection for Ginga is 
given in Hayashida et. al. 1989. 

To help the rejection of the main sources of internal background, the satellite 
recorded some ‘house keeping functions’. Among them: When passing through the 
SAA the detector voltage was lowered, the observation was stopped and the time 
was recorded; Data bins from the next few tens of minutes right after the passage 
through the SAA were discarded because of strong components from the short term 
decay of the activated satellite material. In the other orbits the COR was calculated 
continuously from the local orbit parameters, including the magnetic field, and data 
bins which corresponded with the time in which particles with energy less than 9 
Gev were expected to hit the satellite were discarded; Data bins from parts of the 
orbit in which the source was obscured by earth were discarded, as well as data from 
line of sight which grazed the upper atmosphere and thus contained possible X-ray 
reflection from the sun; The correlation of counts detected in the MID detectors in 
the high energy channels with the internal background (see Hayashida et al. 1989), 
was used for discarding data bins which contained high background; The lower energy 
channels were used for rejection of soft sources, such as random reflections from the 
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upper atmosphere of solar flares. 

For rejection of the external (sky) X-ray background observation of the region 
near the source before and/or after the observation was usually used. For verification 
of this method the method developed by Awaki (1991) is used in our work. 

5.3 Data Analysis 

In order to obtain high statistical significance in the fitting of spectra several 
methods are usually employed. For example, data from the top and mid detector 
layers of LAC are combined for increased detection sensitivity in the 8 — 20keV range. 
Combining these layers results in an increase of some internal background features 
near 3 and 5 keV (which are probably the result of the 41 minutes and 8 hours 
typical decays of satellite materials after the passage through the SAA). For flux 
variability observations eight 16 second collection bins tire combined, and the external 
background is then subtracted. 

Spectral fitting files are accumulated over the whole observation, or in tempo- 
ral groups of certain minimum length so that each spectrum or hardness ratio are 
obtained with high statistical significance. (The hardness ratio, flux in 8 — 20keV di- 
vided by the flux in 2 — 6keV, is an approximate measure of the overall spectral slope.) 
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This technique is used in order to find correlation between spectral features. This 
is a sophisticated temporal analysis of the spectrum which determines the reaction 
times of parts of the spectrum to the overall changes in flux or hardness of flux. The 
changes in hardness of flux are generally used in order to find positive and negative 
delays in the onsetting of changes of the spectral features (Kunieda et al 1990). 

5.4 THE ROSAT MISSION AND PSPC 

ROSAT stands for ROentgen SATellite and is an X-ray observatory developed 
through a cooperative program between Germany, the United States, and the United 
Kingdom. The satellite has been designed and operated by Germany (Trumper 1990). 
It was launched on June 1, 1990. The ROSAT mission began with a six-month, all- 
sky survey, after which the satellite began a series of pointed observations which will 
continue for the duration of the project. 

ROSAT carries three instruments. The main instrumentation consists of a grazing 
angle X-ray telescope (XRT) which can be used with two X-ray detectors: the Posi- 
tion Sensitive Proportional Counter (PSPC) and the High Resolution Imager (HRI). 
The X-ray telescope is complimented by the Wide Field Camera (WFC), an EUV 
telescope. With this instrumentation, ROSAT is capable of observing celestial tar- 
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gets at photon energies in the range from 20 eV to 2 keV (Trumper 1990). Beginning 
in the first half of 1994, all ROSAT-XRT observations were made with the HRI. After 
four years of successful operation, the PSPC instrument ran out of gas and stopped 
operating. The spectral sensitivity of the (late) PSPC is shown in figure 5-3. 

Background rejection for ROSAT is relatively simple when compared with Ginga, 
and is not as critical in the data interpretation stage. We therefore choose not to 
elaborate on background rejection methods for ROSAT (but see Hasinger, Trumper 
and Schmidt 1991 and Nandra and Pounds 1992). The analysis of data from PSPC on 
ROSAT does not require any elaborate methods, such as those mentioned in section 
5.2.3. 
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Chapter 6 

OBSERVATIONS OF 
MCG -6-30-15 

In this chapter we describe the comprehensive analysis of data from the Seyfert galaxy 
MCG -6-30-15. After a brief description of past observations of this and other Seyfert 
1 galaxies we present the observational data from the Ginga and ROSAT satellites. 
-\ We then try to fit the observed flux variability, spectrum and spectral variability with 

current models of AGN central engines. 

6.1 OVERVIEW 

J 

Much has been learned about the nuclei of Seyfert 1 galaxies from their X-ray 
spectra and light curves in the last decade. Ginga and ROSAT missions were instru- 
mental in obtaining that data. Since we used data from these mission we open by 
presenting some of their successful results. 
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The X-ray spectrum of Seyfert 1 nuclei has many features. To first order, the 
2 — 50keV spectrum of the great majority of Seyfert 1 galaxies is dominated by a 
simple power law of the form F(E ) = A.E _r photons cm -2 sec -1 keV -1 (Turner and 
Pounds 1989, Mushotzky, Done and Pounds 1993). This fact, which was known since 
the successful observations by the HE AO-2 satellite (alias Einstein Observatory), was 
verified by Ginga. Observations with Ginga and ROSAT have also revealed additional 
features in the spectra of Seyfert 1 galaxies. Ginga’s broad energy band enabled the 
discovery of a hard ‘tail’ in the 8 — 20keV range (Pounds et al 1990). This tail consists 
of X-ray flux in excess of that expected from the power law emission. The high photon 
count in Ginga was utilized in showing that Fe K-shell emission lines are a common 
feature in the spectra of most Seyfert 1 galaxies (Awaki et al 1991, Nandra et al. 

1991) . Observations with the German satellite ROSAT have revealed photoelectric 
absorption by cold and warm (partially ionized) material in the line sight in the 
0.1 — 2.5keV range. This absorption is larger than that expected from matter in our 
own galaxy, ~ 3 x 10 2 °atoms cm -2 in high galactic latitudes, for most Seyfert Is 
(Pounds et al. 1993). Oxygen absorption edges (for example, see Nandra and Pounds 

1992) , and emission lines (see Netzer, Turner and Pounds 1994 and George, Turner 
and Netzer 1995) were also found in the spectra of Seyferts by ROSAT and the new 


Japanese satellite Asca. 

Most Seyfert 1 nuclei show some evidence of variability. A significant fraction of 
Seyfert Is was shown to have large amplitude ( SF/F > 1) variability on time scales of 
days or less (McHardy 1988, Grandi et al. 1992, Green 1993). Ginga’s good statistics 
helped discover temporal variability on very short timescales, and spectral variability 
on short time scales in most Seyferts. Temporal variability over timescales of hours 
was observed in more than half of all Seyfert Is. Temporal variability over timescales 
of years was found in 97% per cent of all sources. Spectral variability was observed 
in roughly half of the Seyfert Is. The spectrum is variable in the sense that the 
power-law usually steepens when the overall intensity increases (Grandi et al 1992, 
Mushotzky, Done and Pounds 1993). 

We used the advantages of Ginga and ROSAT in our analysis of data from four 
Ginga and one ROSAT observations of MCG -6-30-15. The Seyfert 1 galaxy MCG 
-6-30-15 has a bright X-ray nucleus with relatively high signal to noise ratio. The 
relatively long observations of this source with Ginga made it a good candidate for the 
possible detection with high statistical significance of variability of continuum with 
timescales of down to minutes, variability of the absorption column with timescales 
down to tens of minutes, variability of the hard tail down to hours and variability of 
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the Fe K-shell emission line on time scales down to several hours. 

The results of our analysis of data from MCG -6-30-15 constrain the possible 
models of this source. We verified that the flux from MCG -6-30-15 changed over 
timescales down to tens of minutes, which could be related to a size of the continuum 
emitting region of less than ~ 3 x 10 13 cm. We saw evidence for changes in the 
absorption column over timescales of down to hours. The thickness of this warm 
absorber might therefore be limited to ~ 10 14 cm. The changes in absorption column 
were anti-correlated with the overall continuum flux, and could therefore be attributed 
to warm (partially ionized) matter in the line of sight. We also found evidence that 
the Fe emission line varied on timescales of days, which may correspond with warm 
or cold matter within 10 15 cm of the SMBH. The excess of photons in the hard tail 
was found to be anti- correlated with the continuum flux and correlated with the 
absorption column. There were two possible interpretations of this result. In the first 
an accretion disk Compton scattered continuum radiation that was then absorbed 
by a warm absorber that lied within 10 14 cm from the disk (see Mushotzky, Done 
and Pounds 1993). In the second interpretation the warm absorber itself, which only 
partially covered the continuum source, was responsible for the hard tail. The hard 
tail therefore included photons which were Compton reflected from warm absorber 
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outside the line of sight, and photons which were absorbed by the highly ionized 
warm absorber in the line of sight. In this interpretation the location of the partially 
covering warm absorber could be determined from the variability of the emission lines 
(See Netzer 1993, Sivron et al. 1995). 

Other observations of MCG -6-30-15 agree with the above models. We verified 
the results of ROSAT observations of MCG -6-30-15, which revealed the presence of 
a warm absorber in the line of sight (Nandra and Pounds 1992, Fabian et al. 1994). 
We also used the observation by ROSAT to verify that MCG -6-30-15 field of view 
was free of sources that could have complicated the analysis for the lower resolution 
Ginga detectors. Our results could also be reconciled with recent Asca observations 
(Reynolds and Fabian 1995), but one of our interpretations differs from theirs. 

6.1.1 History of the Observations of MCG -6-30-15 

The existence of cold matter in the central engines of AGN is still an unresolved 
issue. AGN emit strong radiation in wave bands ranging from radio to X-rays. Multi- 
frequency observations have shown that those different wave bands do not originate 
in the same spatial region (Done et. al. 1990). In many of these extra galactic 
sources the spectrum of the X-ray emission itself has complex features superimposed 
upon a power-law feature (Turner and Pounds 1989, Pounds, Nandra and Stewart 
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1992, George and Nandra 1994) which is known to be variable, and thus assumed 
to originate in the central engine of AGN (Guilbert and Rees 1988). Much effort 
has been directed towards finding whether these features also originate in the central 
engine. These features are thought to originate in cold, or partially ionized, matter 
in the central engine. In this chapter we show that in the case of MCG -6-30-15 the 
origin of some of these features is probably very close to the central engine. 

The Seyfert 1 galaxy MCG -6-30-15, at red-shift z = 0.0078, was observed by 
HEAO-2 and EXOSAT before it was observed by Ginga. The detectors on Einstein 
observatory observed MCG -6-30-15 in 1979 and the detectors on Exosat observed 
it in 1984 and 1985. The total luminosity was found to be L x ~ 0.8 x 10 43 and 
Lo. 2-2 ~ 1.2 x 10 43 erg sec -1 , but the luminosity in the EUV range was expected 
to be much larger due to ionization of the Balmer series lines (Reichert et a!. 1985, 
Pounds, Turner and Warwick 1986). An iron emission line was also found by EXOSAT 
(Nandra et al. 1989). Analysis of the temporal behavior of this source by means of 
sophisticated power spectrum techniques have been performed by McHardy in 1988 
and Green in 1993 (see section 1.1.5). These yielded a slope of — 1.35±o25 f° r 
power spectrum for both the Exosat 1985 and the Ginga 1987 observations. 

Ginga observations of MCG -6-30-15 were analyzed in several earlier papers. The 
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Seyfert 1 galactic nucleus of MCG -6-30-15 is among the brightest radio quiet AGN 
observed, which makes the background subtraction easier than in other Ginga weak 
sources (Matsuoka et al. 1990, Nandra, Pounds and Stewart 1990), and the com- 
parison of Ginga observations from four different years therefore reliable. Most of 
the superimposed spectral features mentioned in the previous section are found in 
MCG -6-30-15, which is found to have a variable flux on time scales down to hours 
(Nandra, Pounds and Stewart 1990), and is therefore a prime candidate for the de- 
tection of spectral features variability. Evidences for a ‘soft excess’ below IkeV, a 
strong Fe K - emission line, a possible Fe K - absorption edge and a significant hard 
tail excess in the 10 — 20keV were found in previous analysis of the 1987 GINGA 
observation and the 1984 Exosat observation of this source (Matsuoka et. al. 1990, 
Nandra, Pounds and Stewart 1990). The recent Asca and ROSAT observations have 
revealed the presence of an Oviii absorption edge, and a possible Oyn absorption 
edge (Nandra and Pounds 1992, Fabian et al., 1994, Matsuoka et al. 1994, Reynolds 
and Fabian 1995). Data from the ROSAT and Asca X-ray telescopes have also shown 
that the absorption by cold matter in the line of sight could not be much more than 
the galactic value. The recent observation of this source with Asca in the 0.5 — lOkeV 
energy band, a little bit lower than that of Ginga, has shown it to be spatially clean of 
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interferences from other sources (especially ones with similar spectrum) in a 1 f.o.v. 
(Nandra and Pounds 1992, Fabian et al. 1994, Inoue 1994), which makes it free of 
possible confusion with galactic sources. 

Our analysis of MCG -6-30-15 focuses on problems that were not addressed in 
previous papers. Although the flux variability and spectral variability of this source 
were analyzed by previous authors (September 1987 data, Matsuoka et. al. 1990, 
Nandra Pounds and Stewart 1990, Fiore et al. 1992) we note that their analyses 
did not take advantage of the low noise in the late observations data (June 1989, 
February 1990, July 1990) by GINGA, which we extensively use in this thesis. The 
transmission rate in those later observations is much improved when compared with 
the 1987 observation. The main focus of this chapter is much different from the earlier 
papers in that the emphasis is on variability of the spectral components over ranges 
from minutes to years, and thus a maximum utilization of the multiple observations 
is achieved. The fact that the source was observed with the same instrument is 
extremely helpful in establishing year to year variations. 

6.2 OBSERVATIONS 

In this section we describe the observations of MCG -6-30-15 by the LAC detectors 


on Ginga and the PSPC detectors on ROSAT. The LAC and the PSPC are described 
in chapter 5. 

The conditions of the LAC observations were determined by machine and orbit 
constraints. The LAC observations were performed in the MPC-1 mode, in which 
the X-ray signals were accumulated in 48 energy-sorted channels. The temporal 
resolution in this mode is 16 seconds. To obtain better statistical significance in 
both the temporal and spectral studies of this source we combined the top and mid 
layers. This method was effective for the later observations, in which the internal 
background was reduced due to the decay of Ginga’s orbit. The accumulation time 
for the light-curves was chosen to be 128 seconds, just above the Poissonian noise. 
The maximum continual observation time obtainable by Ginga is roughly 1000 — 2000 
seconds. Interruptions due to Earth occultation, passage through the South Atlantic 
Anomaly (SAA) and other high background regions can be noticed in the light curve 
in Figure 6-2. A 1000 seconds accumulation time is also the shortest for a spectrum 
of this source. The longest continuous observation time we used was five, almost 
consecutive, observation days (see table 6-1). 



Table 2. MCG -6-30-15 observations in this thesis: 


File names'* 

Date of Observation 1 * 

Instrument' 

Flux d 

87 

October 9-10 1987 

LAC 

io- u 

89 

June 25-26 1989 

LAC 

10" n 

89a° 

89lw,md,hi 

891mdmin,himin <s 

89sft,hrd 

90 

February 2-7 1990 

LAC 


90a e 

90lw,md,hi 

901wmin,mdmin,himin 

90sft,hrd 

90dl 

February 2 



90d2 

February 3 



90d3 

February 5 



90d4 

February 6 



90d5 

February 7 




a File names after comma have the same first three letters. 
b The dates in which the observation took place. 

c LAC: Large Area Counter on Ginga. PSPC: Position Sensitive 


Proportional Counter on Rosat. 
d Flux at 1 Kev. 



138 


O 


O 


r 


c 


e 


o 


c 


G 


O 


o 


MCG -6-30-15 (a = 13 A 33 m 00 4 , S = — 34 o 02 , ) was observed on four separate occa- 
sions with Ginga. In the September 10, 1987 observation two of eight detectors, placed 
at an end of the 4x2 detector array, were excluded because they were illuminated by 
the sun at angles < 75°. The other detectors were in the shadow of the solar panel. 
Only top layer data was used, because at that early stage of the mission the internal 
background was the largest, and because the mid layer exhibits a lower efficiency for 
X-rays when compared with the particle noise. Also, the collimator efficiency was 
mediocre (< 85%, Nandra, Pounds and Stewart 1990). These handicaps resulted in 
different values for the same spectral features when different background substruc- 
tion were used (See Nandra et.al. 1990 for details). In reanalyzing this measurement 
we add the mid layer, which significantly increases the efficiency in the 10 — 20keV 
range (see section §4). Adding the midlayer increases the internal background be- 
cause the satellite was at a high altitude in the early stage of the mission, but must 
be done in order to be consistent with the later observations. The background sub- 
traction method, similar to that performed by Matsuoka et.al. 1990, used nearby sky 
background, observed within several days, in the direction a = 13 fc 13 m , — 38°42\ 

In the June 25-June 26, 1989 observation the collimator efficiency was high, the 
overall background was lower because of the satellite descent into the upper atmo- 
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sphere and the low solar activity, and the adjacent day background observation was 
in MPC-1 mode. The overall flux level was higher and all detectors were pointed at 
an angle > 90° from the sun. We therefore used an off source data base made from 
observations in the direction (a = 8 = — 35°58 ) on June 24, 1989. 

By far the longest observation of MCG -6-30-15 with Ginga occurred on February 
2nd through 7th 1990. During that period the LAC field of view drifted on several 
occasions, resulting in low collimator efficiency during most of the February 4th obser- 
vation, and the last few orbits on February 2nd and 5th, all of which were discarded. 
For background subtraction two methods were used. A nearby observation of blank 
sky on February 12, 1990, in the direction a = 13*31”*, S = -35°58' was used as 
background for the observations. We also used the sky averaging method described 
by Awaki 1991. The nearby blank sky observation was chosen, because of better cor- 
relation of mid to top layer ratio with the theoretical value. In figure 6-1 we see the 
correlation of the top to mid ratio with its theoretical value for the 1990 observation. 


EE 


energy c KEV j 


Figure 6.1: The expected theoretical ratio of the counts observed in 
the top layer of the LAC detectors over the mid layer is shown in solid 
steplike curve. The observed value is shown by the error bars. This 
method is used for verification of data quality. 
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We have reanalized the PSPC observation of MCG -6-30-15 on January 29, 1992 
using standard background subtraction and correcting for the 10 arc-minutes mis- 
alignment which can be seen in figure 6-2. In figure 6-2 we see that MCG -6-30-15 
was roughly 10 arc-minutes off target on the date of the observation. It is also clear 
that there are no major interferences from other sources in that observation. 
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Figure 6.2: Rosat observation of MCG -6-30-15 on January 1992 shows 
the field of view to be free of intense hard sources. 
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6.3 INTENSITY VARIATIONS 

Before looking at the X-ray spectrum of an AGN observers usually consider the light 
curve in the different wavebands. Visual inspection of these light curves usually 
indicates the general spectral variability behavior of the source. 

6.3.1 Light Curves 

The light curves observed by the top and mid layers of LAC in the first day of the 
1990 observation are presented in figure 6-3. The total (background subtracted) flux 
of 2 - 20keV photons is shown in the lowest panel. The lower energy cutoff (2keV) 
of the low band, shown in the middle panel, is chosen so that the soft excess, and 
the Ovm edge of MCG -6-30-15 (Nandra and Pounds 1992, Walter and Fink 1993) 
are not included. The high energy cutoff (4keV) of the low energy band was chosen 
so that the Fe K line (6.4keV) is not be included. The low energy cutoff (8keV) of 
the high energy band was chosen to avoid possible contribution from the Fe edge 
(7.1 - 8.8keV, depending on ionization state, but known to be less than 8keV in this 
source), if such edge existed. The high energy cutoff (20keV) of the high energy band, 
shown in the upper panel, was chosen to discard the energy range with low signal to 


noise ratio. 


144 


The light curves in the 1987, 1989 and 1991 observations show similar behavior 


(however, see chapter 7). 


<*«V /l-ZohV 



Figure 6.2: The flux variability observed by Ginga in the first two days of the 1990 obser- 
vation. The variability of the total flux is shown in the lower panel, the variability of the 
2 — 4Kev flux is shown in the second panel. The variability of ho 8 20Kev flux is shown 
in the third panel. The ratio of the 8 - 20Kev over the 2 - 4Xev flux is shown in the top 
panel. The abscissa unit ij seconds. 
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6.3.2 Correlation of the Two Energy Bands 

In the upper panel of Figure 6-3 we see correlation between the low and high 
energy bands. The hardening of the spectrum with decreasing flux is obviously present 
throughout the observation period. No delays can be inspected by eye (delays with low 
significance level were found in Piro et al. 1992 for the 1987 observation). Assuming 
that specific spectral features that were emitted from the AGN are responsible for 
these changes we can differentiate between models in which the hardening with flux 
is expected and other models. 

In figure 6-4 we see the correlation of the 8 — 20keV and 2 — 4keV energy bands 
with 128 seconds time bins. Note that the data cannot be fit with a line that goes 
through the origin. The existence of hard photons in the zero soft photons limit may 
be attributed to reflection or partial covering of a power-law spectrum (see section 
6.4), whereas the scattering may be attributed to delays. 
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6.4 CHANGES IN SPECTRAL FEATURES 

The various parts of the spectrum may be related to processes in the nuclei of MCG 
-6-30-15. In this section we present and evaluate models that may explain the spectral 
features. 

6.4.1 Preliminary inspection of the sorted spectrum 

In figures 6-5 and 6-6 we see the raw data after background subtraction for the 1990 
observation. In figure 6-5 the accumulation time is the entire observation, whereas in 
figure 6-6 the data was accumulated by flux level with binning time of roughly tens 
of minutes. The data was accumulated by flux level with binning times of minutes, 
tens of minutes, days and years (87, 89, 90 and 91, see table 2). An accumulation of 
data by hardness ratio with binning time of tens of minutes, as described in chapter 
5, results in spectra which are similar to the low and middle flux levels in figure 6-6. 
Figure 6-7 includes the 1992 observation by PSPC on ROSAT. Similar spectral files 
for the 1987, 1989 and 1991 observations are described in Sivron et al. 1995. 
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Although some of the decrease in photon count in the ~ 2 — 3keV range can be 
attributed to detector efficiency (see figures 5-2 and 5-3 in chapter 5), it is evident 
that some of the decrease is due to true observed spectrum. This decrease is clearly 
anti correlated with flux, and also anti correlated with softness. In the 3 — 5.5keV 
the photon counts spectrum seems to have a negative slope, which decreases in the 
8— 14keV range. The decrease is anti correlated with flux and softness. In the ROSAT 
raw data there is deficiency in counts in the 0.2 — 0.9keV range, as compared with 
the detector efficiency and power-law emission (see section 6.4.6 below). 

6.4.2 The Absorption Column 

The absorption column of matter in the line of sight to MCG -6-30-15 is probably 
smaller than the Nh ~ 10 21 ' 6 cm -2 value previously considered for the Ginga obser- 
vations (Nandra, Pounds and Stewart 1990). Preliminary results from ROSAT and 
ASCA find N h 10 20 7 cm -2 with high statistical significance (Nandra 1991, Fabian 
et. al. 1994). Earlier observations with MPC on Einstein observatory (Reichard 
et.al.1985) and ME combined with LE on EXOSAT (Nandra et al. 1989) are consis- 
tent with the smaller values. The approximate galactic value for the cold absorption 
is roughly Nh = 10 2a7 cm -2 (Walter and Fink 1993 and references therein). 

It may be noted that fitting models which include warm absorption in or near 


the central engine may introduce an additional absorption column as an artifact of 
the multi-parameter fit, especially as a result of bound-free transitions in partially 
ionized oxygen. We therefore set the value of the absorption column at its galactic 
value, for most of our spectral fits. In figure 6-8 a schematic picture of the effect of 
absorbing material in the line of sight on a power-law spectrum is shown. 
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Figure 6.8: A schematic picture of the effect of absorbing material in 
the line of sight on a power-law spectrum. The upper curve in the 
graph is the unabsorbed spectrum. 



6.4.3 Power law with an absorption feature and an Fe emis- 
sion line 

The simplest description of the above spectra is with 

F(E) = (N x E~ r + iV Fe e" ((E - E# '‘ )J/<7 ^ ) )e ( - <,(E)xWw) , (6.1) 

where N and Np e are the normalizations for the power law and Gaussian emission 
line, E is the photon energy, a(E ) is the cross-section for cold cosmic abundance 
matter, Ep e is the iron emission line energy and <r Fe its observed width. 

We fitted the spectra with an iron emission line with local frame energy 6.4keV, a 
power-law spectrum with photon index in the 1.8 — 2.2 range, and a cold absorption 
column in the range Njj — 10 2 ^~ 22 cm *. The natural width of the line was assumed 
to be 0.05keV, but after Fabian et.al. 1994 reported an intrinsic width of ~ 0.4keV 
the data was fit again, albeit with negligible changes, because the energy resolution 
of the LAC proportional counters is IkeV. 

The results from the 1990 observations are summarized in tables 3, 4, 5 and 6. 
The fit is unacceptable in most cases, and we the residuals, using the 1990 file in 
figure 6-9, show an edge-like feature near 8keV with excess of hard X-ray counts, 
usually referred to as the hard tail. 

Some results can be derived by fitting together spectral accumulation files for 
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different flux levels in different energy bins. In tables 3, 4 and 5 one can see that 
the normalization of the power-law emission is variable with at least 99.99% confi- 
dence level (see Bevington 1968). The absorption column is also variable with 99.9% 

confidence level. We assumed this variable absorption column to be an indication 

C 

of internal changes inside the source, which can probably be associated with vari- 
able warm absorption or variable covering factor (see section 6.4.6). In figure 6-10 
C we see that the photon index T and the absorption column are strongly correlated, 

suggesting that these parameters cannot be individually determined with the above 
model. 

It is not at all clear that the Fe line flux is variable for all time binnings. The 
confidence level is best for binning by days (90 %) and binning by years (95 %), 
both of which are only marginally acceptable. Figure 6-11 is a schematic picture that 
describes the possible location of the Fe line emitting matter. 
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Figure 6./$ The x contour plot for the absorption column and the 
photon index in the second day of the 1990 observation. The inner 
contour corresponds to 99% confidence level of fit, the second contour 
corresponds to 95 % confidence level, and the outer contour corresponds 
to 67 % confidence level. 
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Figure 6.11: A schematic picture that describes the possible location 
of the Fe line emitting matter around the power law source. The graph 
describes the subsequent delay in Fe emition flux with respect to the 
power law flux. 





6,4.4 The Hard tail and the Apparent hard ‘Edge’ 


There are two models with which we tried to determine if the hard tail and edge 
features are variable. The first one is the disk reflection model. The disk model can 
be written as 


F(E) = (N x £' r x [1 + rA(E)] + N Fc e tiE - B ” , ' l ’r- ) )e l 


(6.2) 


with r = fi/ 47 r, where is the solid angle subtended by the disk, and A(E) is the 
Albedo as calculated, for example, in appendix . 

The second model is the partial covering by cold absorber model, or the leaky 
absorber model. The leaky absorber model can be written as: 


F(E) = (N x £- r x [1 + f c J{ 1 - f cov )e~°W N "‘~] + 


(6.3) 


where f cov is the covering factor, and N Hc ov is the absorption column in the cold matter 
which partially covers the source. We note that for column densities in excess of 10 23 
atoms cm -2 this model is not accurate in a three dimensional manner, because some 
reflection is expected from the absorbers on the other side of the power-law source. 
For covering factor of much less than ~ 0.5 this contribution is probably very small. 
(See also discussion in chapter 7.) 
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The improvement in the fits is similar with both models (see tables 3 through 6). 
There is also a similar improvement in confidence level when the solid angle or the 
covering factor become free components. When we fit the different flux level files and 
let / vary we get significant improvement in fits of about 99.99% significance level in 
the days and years binning. For r we get 99.99% significance level of improvement for 
all binnings. This is incredibly important, because it implies that the disk changes 
on time scales of minutes. In order to verify this point we also checked the files 90hrd 
and 90sft, with similar results. 

Schematic pictures of the reflection and leaky absorber model axe shown in figures 
6-12 and 6-14. The fits and residuals in both of the above cases are shown in figures 
6-13 and 6-15. 
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Figure 6.14: A schematic picture of the leaky absorber model, and the 
resulting spectrum. 
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6.4.5 Soft X-ray features 

Since the soft x-ray absorption like feature is correlated with the hard tail we 
considered models that naturally relate these two. The first one is the disk plus 
warm absorber model. To show that warm absorption is indeed possible we fit the 
ROSAT 0.1 — 2.5keV data with an edge feature of the sort 

F(E) = N x E~ r x e -HE-E{ e d ge )) x e [-<r(E)xN„)^ (6 .4) 

where a(E) is set to the galactic value, and t(E - E(edge )) is the edge optical depth. 
In figure 6-16 an edge is indeed apparent in rest frame energy 7.8keV, corresponding 
to highly ionized Ovii or Oviii oxygen (see Nandra 1992). 
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We tried to tie together the coefficients responsible for soft absorption and hard 
tail in the Ginga data. The fit improves with high confidence level. In order to include 
more general cases, in which the warm absorber may be, for example, highly ionized 
in one distance and less ionized in another, we modeled the disk reflection as we did 
before, but with variable edge (see table 6). In figure 6-16 a schematic picture of the 
disk plus warm absorber model is presented. 

The partial covering model of the previous section was approximated as 

F(E) = (Wx£- r x[l+/ lcov /(l-/icov)e _<r(E)JVwtowl ] 

(6.5) 

x [1 + f 2eov /( 1 - h cov )e- a(E)NlUov '] + N Fe e^ E ~ E ^ 2 ^), 

for small covering factors. There is hardly a case in which there is significant improve- 
ment when the first covering factor becomes different from the second covering factor, 
but the mere addition of another absorption column improves the fits tremendously 
in all time bins (see tables 3, 4 and 5). In figure 6-17 a schematic picture of the partial 
covering by a warm absorber is presented. 
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Table 3. Hours, Partial Covering,90 


90LW 1.5 ±0.09 1.82 ±0.03 0.34 ±0.12 1.04 ±0.49 
2.5 1.82 ± 0.03 0.34 ±0.12 1.6 ±3.0 

2.6 ±0.8 1.93 ±0.10 4.3 ±2.9 1.4 ±0.3 


1.4 ±0.3 0.29 ± 0.09 340 ± 200 0.97/0.53 


6.0 ±7.3 

2.3 ± 0.4 

3.1 ±2.2 

0.5 

0.58/0.48 

54 ±20 

0.89/3.6 

2.6 ±0.07 

1.92 ±0.01 

0.34 ± 0.05 

1.4 ±0.3 

... 

... 

1.29 

4.1 

1.92 

0.34 ± 0.05 

2.2 

0.38 

1000 

1.4/0.01 

4.7 ±0.6 

2.03 ±0.05 

3.4 ± 1.1 

2.2 ± 0.5 

0.33±).04 

310 ±80 

1.05/8.2 

4.5 ± 1.0 

2.05 ± 0.08 

3.0 ±1.1 

2.0 ± 0.5 

0.38/0.26 

190 ± 100 

0.96/3.3 

3.3 ±0.1 

1.94 ±0.02 

0.35 ± 0.08 

1.2 ±0.6 

... 


1.07 

7.4 

1.96 ±0.05 

0.36 ±0.13 

2.7 

0.55 

1000 

1.15/0.14 

6.2 

1.96 

0.99 

2.2 

0.45 ± 0.5 

1000 

1.16/0 

10.9 ± 10 

2.35 ± 0.33 

3.15 ± 1.44 

0.9 

0.55/0.42 

70 ±10 

1.14/1.3 


1.4 ±0.5 
1.4 ±0.3 
1.0 ± 0.6 
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Table 3 — Continued 


File 

Norm* 

r 


/f. c 

f d 

Nh7° 

X%/F ( 


6.6 ±3.1 

2.02 ± 0.4 

1.2 ±1.2 

2.9 ± 2.5 

0.49 ± 0.27 

1000 

1.63/0.0 


11.0 ±8.9 

2.44 ± 0.29 

2.5 ±1.2 

1.0 ±1.2 

0.64/0.44 

55 ±13 

1.47/3.2 

90HRD 

1.9 ±0.1 

1.81 ±0.03 

0.27 ± 0.08 

1.5 ±0.4 

... 

... 

1.33 


3.1 ± 1.5 

1.81 ±0.03 

0.27 ±0.08 

2.5 

0.39 

1000 

1.44/0.0 


3.1 ±0.7 

1.89 ±0.09 

2.9 ±2.1 

2.2 ±0.7 

0.28 ± 0.07 

260 ± 160 

1.39/0.9 


5.8 ±5.9 

2.19 ±0.36 

3.3 ± 1.5 

1.6 ±1.0 

0.50/0.40 

70 ±17 

1.39/1.1 


varpo 2.92 ± 0.09 
2.28 ±0.07 


2.72 


varFe 


1.00 ± 0.45 

2.13 ±0.39 2.57/4.24 


varf 


2.00 

1.84 


1.44/44.7 


varfl 0.0 

i 

0.27 ±0.03 1.44/46.9 


vHl 2.02 


1.95 


1.45/0.72 
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Table 3 — Continued 


File 

Norm 4 T Nn h If* c 

/ d 

NhS xllF 

flf2 

set 

0.58 



set 

0.36 

1.49/0.7 

vflf2 


0.64/0.50 

0.68/0.62 

1. 4/3.8 


*in lG" 2 erg sec 1 cm" 2 at lKev. 
b Cold Absorptions, 10 22 cm" 1 
c in 10~ 2 erg sec 1 cm" 2 at lKev. 

^Covering factor/second covering factor. 
e Cold Absorptions, lO^cm" 1 , coefficient 2. 

f Note: 31 channels per file. The F test: level of confidence in an additional parameter. 
5 Power law normalization is variable in different files in the fitting program 
h Variable Gaussian norm 
‘Variable Photon index 

^Variable one partially covering cold absorber covering factor 
k PC Cold absorber + variable photon index 
l Two PC cold absorbers, with a single covering factor 
m Two PC with two covering factors 


References. — (1) me. 
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Table 4. Days, Partial Covering, 90 


File 

Norm* 

r 

cr 

h.' 

f d 

N„,' 

X 2 r /F ( 

90D2 

2.1 ±0.1 

1.95 ±0.03 

0.57 ±0.10 

1.7 ± 0.4 

. . . 


0.82 


4.6 

1.95 ±0.09 

0.59 ±0.19 

3.5 

0.5 

1000 

0.38/0.17 

... 

4.6 ± 1.5 

1.98 ±0.04 

1.6 ± 1.4 

3.3 ±2.1 

0.51 ± 0.20 

1000 

0.86/0.6 

... 

3.2 ±1.0 

2.04 ±0.13 

2.4 ±1.7 

2.1 ± 0.9 

0.48/0.21 

220 ± 290 

0.83/2.2 

90D3 

2.5 ± 0,1 

1.88 ±0.02 

0.23 ± 0.07 

1.3 ±0.4 

... 

... 

1.04 

... 

3.9 ±0.9 

1.92 ±0.03 

0.33 ± 0.08 

2.1 ± 1.0 

0.33 ±0.16 

460 ± 320 

0.83/8.8 

... 

4.5 ±1.0 

2.02 ± 0.07 

4.2 ±1.7 

1.3 ±0.7 

0.29 ± 0.06 

260 ± 100 

0.72/3.76 

... 

4.3 ±1.0 

2.0 ± 0.06 

4.3 ± 1.8 

2.0 ± 0.8 

0.27/0.35 

350 ± 180 

0.70/0.62 

90D4 

3.0 ±0.1 

1.88 ±0.02 

0.16 ±0.08' 

1.3 ± 0.6 

... 

... 

1.67 

... 

5.4 ±1.2 

1.94 ±0.03 

0.28 ±0.09 

2.6 ± 1.3 

0.40 ± 0.14 

430 ± 230 

1.37/8.07 

... 

4.7 ± 1.1 

1.94 ± 0.06“ 

0.73 ±0.16 

3.7 ±1.7 

0.5 ±0.1 

330 ± 160 

1.36/0.07 

... 

5.7 ±1.3 

1.96 ± 0.06 

1.3 ±2.4 

2.6 ± 1.2 

0.31/0.41 

420 ± 230 

1.42/0.01 

90D5 

2.5 ±0.1 

1.88 ±0.02 

0.23 ± 0.07 

1.3 ± 0.4 

... 

... 

1.04 . 


3.9 ±0.9 

1.92 ±0.03 

0.33 ± 0.08 

2.1 ± 1.0 

0.33 ±0.16 

460 ± 320 

0.83/8.8 


4.5 ± 1.0 

2.02 ± 0.07 

4.2 ± 1.7 

1.3 ±0.7 

0.29 ± 0.06 

260 ± 100 

0.72/3.75- 

... 

4.8 ±1.0 

2.0 ±0.06 

4.3 ±1.8 

2.0 ± 0.8 

0.27/0.35 

350 ± 180 

0.70/0.62 

90D6 

3.0 ±0.1 

1.88 ±0.02 

0.16 ±0.08 

1.3 ± 0.6 

... 


1.67 

... 

5.4 ±1.2 

1.94 ±0.03 

0.28 ± 0.09 

2.6 ± 1.3 

0.40 ±0.14 

430 ± 230 

1.37/8.07 

... 

4.7 ± 1.1 

1.94 ±0.06 

0.73 ±0.16 

3.7 ±1.7 

0.5 ±0.1 

330 ± 160 

1.36/0.07 


5.7 ±1.3 

1.96 ±0.06 

1.3 ±2.4 

2.6 ±1.2 

0.31/0.41 

420 ± 230 

1.42/0.01 
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1.96 


1.07/35.4. 








2.5 ± 0.1 1.34 ± 0.02 

4.4 ± 1.3 1.91 ± 0.02 

4.1 ± 0.7 1.^3 ± 0.4 

4.7 ± 1.1 1.96 ± 0.05 

1.7 ± 0.1 1.77 ± 0.04 

4.0 ± 1.0 1.91 ± 0.09 

12.0 ± 7.0 2.3 ± 0.2 

12.1 ± 7.4 2.3 ± 0.2 


0.20 ± 0.06 1.3 ± 0.4 
0.29 ± 0.07 2.4 ± 1.3 

1.4 ± 1.0 2.1 ± 0.9 
3.2 ± 1.9 2.1 ± 1.0 
0.50 ± 0.16 1.4 ± 0.7 
0.82 ± 0.23 3.3 ± 1.7 

6.6 ± 1.6 4.1 ± 2.1 

6.6 ± 1.6 4.0 ± 2.2 


0.41 ± 0.18 610 ±340 
0.33 ± 0.11 420 ±250 
0 . 25 / 0.38 450 ± 250 

0.46 ± 0.11 270 ±150 

0.61 ± 0.09 200 ± 50 

0 . 62 / 0.61 200 ± 70 


1.14 

0 . 84 / 11.5 

0 . 31 / 0.95 

0 . 31 / 0.93 

1.45 

1 . 16 / 8.8 

0 . 30 / 11.3 

0 . 83 / 0.01 


2.4 ± 0.07 
1.9 ± 0.06 







a 
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1.3 ± 0.1 1.73 ± 0.04 1.1 ± 0.18 1.04 ± 0.56 ••• ••• ^ 

14.0 1.85 ± 0.11 1.4 ± 0.3 1.04 ± 1.9 0.89 ± 0.21 1000 1 . 63 / 8.8 

10.5 ± 11.1 1.35 ± 0.05 1.8 ± 0.6 7.5 ± 9.3 0.85 ± 0.16 1000 1 . 61 / 0.3 

20.0 ± 30.0 1.94 ± 0.10 4.2 ± 1.4 10 ± 13 0 . 64 / 0.90 1000 1 . 6 / 1.2 

2.66 ± 0.08 1.91 ± 0.02 0.24 ± 0.06 1.3 ± 0.4 . ••• 1-18 

3.9 ± 0.5 1.96 ± 0.03 0.35 ± 0.08 2.0 ± 0.7 0.27 ± 0.09 310 ± 180 0 . 88 / 11.3 

4.5 ± 0.7 2.00 ± 0.06 2.3 ± 1.4 2.1 ± 0.7 0 : 31 ± 0.05 290 ± 120 0 . 86 / 0.73 

4.3 ± 0.7 2.01 ± 0.06 1.9 ± 1.5 2.0 ± 0.7 0 . 34 / 0.28 240 ± 140 0 . 88 / 0.0 

3.0 ± 0.1 1.88 ± 0.02 0.16 ± 0.08 1.3 ± 0.6 1-67 

5.4 ± 1.2 1.94 ± 0.03 0.28 ± 0.09 2.6 ± 1 . 3 . 0.40 ± 0.14 430 ±230 1 . 37 / 8 . 0 / 

4.7 ± 1.1 1.94 ± 0.06 0.73 ± 0.16 3.7 ± 1.7 0.5 ± 0.1 330 ± 160 1 . 36 / 0.07 

5.7 ± 1.3 1.96 ± 0.06 1.3 ± 2.4 2.6 ± 1.2 0 . 31 / 0.41 420 ±230 1 . 42 / 0.01 

2.5 ± 0.1 1.88 ± 0.02 0.23 ± 0.07 1.3 ± 0.4 ••• 104 

3.9 ± 0.9 1.92 ± 0.03 0.33 ± 0.08 2.1 ± 1.0 0.33 ± 0.16 460 ±320 0 . 83 / 8.8 

4.5 ± 1.0 2.02 ± 0.07 4.2 ± 1.7 1.3 ± 0.7 0.29 ± 0.06 260 ± 100 0 . 72 / 3.76 

4.8 ± 1.0 2.0 ± 0.06 4.3 ± 1.8 2.0 ± 0.8 0 . 27 / 0.35 350 ±180 0 . 70 / 0.62 

3.0 ± 0.1 1.88 ± 0.02 0.16 ± 0.08 1.3 ± 0.6 ••• 1-67 

5.4 ± 1.2 1.94 ± 0.03 0.28 ± 0.09 2.6 ± 1.3 0.40 ± 0.14 430 ±230 1 . 37 / 8.07 

4.7 ± 1.1 1.94 ± 0.06 0.73 ± 0.16 3.7 ± 1.7 0.5 ± 0.1 330 ± 160 . 1 . 36 / 0.07 

5.7 ± 1.3 1.96 ± 0.06 1.3 ± 2.4 2.6 ± 1.2 0 . 31 / 0.41 420 ± 230 1 . 42 / 0.01 


varpo 1.98 ± 0.05 



Table 5 — Continued 


File Norm* V Nn b If* / d NW 


2.42 ± 0.06 
2.93 ± 0.07 
2.42 ± 0.06 
2.93 ± 0.07 

varFe 


varr 1.76 

1.88 

1.76 

1.88 

1.9 

varfl 


2.1 ± 0.6 
1.8 ± 0.4 
2.1 ±0.6 
1.8 ±0.4 
1.0 ±0.6 


0.57 ± 0.07 
0.28 ± 0.07 
0.57 ± 0.07 
0.28 ± 0.07 


XU? 


1.78/3.64 


1.41/24.6 


0.3 ±0.08 


14.3 


c 
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Table 5— Continued 


File 

Norm* T Ifc c 

f* Ntn' 

xyp 



1.34 





1.93 





1.96 


1.07/35.4 


flf2 

set 

0.49 




set 

0.31 




set 

0.49 




set 

0.31 




set 

0.28 

0.92/10.9 


vflf2 


0.43/0.54 





0.29/0.36 





0.43/0.54 





0.29/0.36 





0.25/0.4 

0.9/4.2 



•See table notes, table 3 
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Table 6. Hours, plref 


File 

Norm* 

r b 

Nh' 

If. a 

n/4«- 

X 2 h/P 

S90LW 

1.6 ±0.2 

1.9 ±0.1 

0.45 ±0.19 

1.0 ±0.5 

0.5 ±0.7 

0.93 

S90MD 

3.0 ±0.2 

2.04 ± 0.48 

0.52 ± 0.09 

1.4 ±0.3 

0.9 ± 0.4 

1.03 

S90HI 

3.6 ±0.3 

2.00 ± 0.07 

0.43 ±0.13 

1.2 ± 0.6 

0.4 ±0.5 

1.87/24 

varpo* 

1.9 ±0.1 







2.9 ±0.1 







3.6 ±0.2 





1.29 


‘See table notes, table 3 
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6.5 DISCUSSION 

Important Questions and answers about variable features in MCG -6-30-15: 

1. Is the basic power-law norm variable? The answer is yes for all of the models 
applied in this thesis. Is the photon index variable? In the partial covering by a warm 
absorber the answer is probably not, because not only is there more improvement for 
two absorption columns but also the changes in the photon index are small (< 0.15). 
For the disk the changes in the power law index are also very small, and variations 
of fl are much more probable. 

2. Do the superimposed absorption/reflection spectral features vary? They certainly 
do! For the two absorption column case the variability of the higher absorption 
column is questionable on very short time scales, but certain on time scales of hours, 
days and years. 

2. In the case of variable partial covering: Can we find any reliable relation between 
absorption and partial covering fraction? A positive answer may strengthen the case 
for partial covering vs. warm absorber (see also appendix F). 


Chapter 7 

DISCUSSION AND 
CONCLUDING REMARKS 


In this chapter we present some concluding remarks on this thesis. In the first section 
the two component cloud model is considered, and compared with the more ‘standard’ 
disk-corona model. In the second section we discuss the dynamics of the central engine 
and the merits of the shock model. In the last section we discuss the results of our 
analysis of data from MCG -6-30-15. 

7.1 THE CLOUD MODEL 

7.1.1 Accretion Disks Versus Clouds 

In many cases the two component cloud model is advantageous over the disk- 
corona model. We argue that, at least in the case of high accretion rates, cloud 
configuration is more plausible than disk configuration. We first present some theo- 
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retical arguments, followed by evidence from observations. 

The main advantage of an accretion disk is its ability to transfer angular momen- 
tum outwards and its axisymmetry. This ability is of less important in the central 
engine, where the typical densities and magnetic fields are sufficient for the outwards 
transfer of angular momentum in semi- spherical configurations. 

One theoretical reason to consider inhomogeneities, such as clouds, is the fact that 
the standard thin disks assume that the shear and the resultant angular momentum 
transfer happens in convection cells of size which is proportional to the height of the 
disk. In the central engine such assumption means that the cells are comparable in 
size to the size of the source itself! This picture is more apparent in chapter 4, in 
which an optically thin disk is perturbed by an inhomogeneity. In this thesis we also 
note that all of models of variability require inhomogeneities and strong magnetic 
fields (see chapters 2 and 4). 

Significant emission in the UV and EUV parts of the spectrum may be easier to 
explain with clouds. It was shown that UV emitting accretion disks may be unstable 
in the central engine. With this result in mind it is hard to explain the UV band, 
which contains a significant portion of the total emitted radiation, which is supposed 
to have originated in the central engine. In order to produce so much UV a large 
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surface has to absorb and re-emit the X-rays and 7-rays from the central engine (see 
chapters 2, 3 and 4). Very small amount of matter in the form of clouds is needed 
in order to reprocess the radiation from the central engine. In order to reprocess a 
significant portion of the emission from the central engine an increased surface area is 
needed (see chapter 3). If the surface is not large the reprocessing probably results in 
the same thermal instability which disrupts the disk, especially when a small portion 
of the disk is illuminated by a flux doubling flare. 

Ultimately the alternative geometrical configurations may be distinguished with 
X-ray telescopes with very high resolution. Present day telescopes, however, cannot 
achieve this resolution. In the outer regions of AGN the resolution of the NLR and the 
molecular torus is within reach with optical telescopes of sub-arcsecond resolution. 
The HST observation of a non-axisymmetric accretion disk in a radio galaxy already 
shed some light on the configuration. (The observed configuration already posses some 
problems for the standard accretion disk, because viscous instabilities are expected 
with varying M. See, for example, Chakrabarti 1995.) But the resolution is now 
limited to arc-minutes in the UV to X-ray regime (ROSAT HRI), in which most of 
the energy is emitted. Improvement of at least two orders of magnitude in resolution 
without compromising the spectral resolution is thus needed in order to resolve the 
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geometrical configuration of the central regions. Until such technology is achieved the 
best way to determine the geometrical configuration is through temporal and spectral 
analysis of data. The recent campaigns of multi wave-band observations by Tsuruta 
et al. and others may thus prove to be invaluable in distinguishing between existing 
models. 

One way by which we can distinguish between the disk and cloud models come 
from the observation of polarization. In thin disks scattering atmospheres should 
induce a high degree of polarization (Netzer 1990, and references therein). On the 
other hand, our cloud model predicts very low polarization, because for the size 
of clouds of our kind the reflection directionality in Equation (3.10) results in weak 
polarization (Sivron and Tsuruta 1993). Observations favor weak polarization (Netzer 
1990). Polarization in the X-ray regime, which can improve our determination of the 
geometry, may be achieved in the near future. 

No sharp Lyman edges (as predicted by many disk models) are expected to be 
observable with the cloud model, because the strong feature emitted by the clouds 
is broadened and smeared by e.g. the Doppler shifts with a variety of cloud veloc- 
ities. This is consistent with observations (Antonucci, Kinney, and Ford 1989, and 
references therein). 
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The spectral behavior of the NGC 4051 data near the B-band was examined (Done 
et al. 1990) which shows no rapid variability during the X-ray flares, and find that the 
low energy end of the UV bump already appears there. Therefore the lack of rapid 
variability at this wavelength is consistent with our cloud model. In our model the 
power-law portion of the IR-optical radiation is due to^different mechanisms (Tsuruta 
1990), and hence no correlation is predicted between the IR-optical-UV bands and 
the X-ray band. On the other hand definite correlated variations are predicted among 
the IR-optical-UV and X-ray bands when a variety of disk models are adopted; for 
instance, the standard thin-disk model (Shakura and Sunyaev 1973) and a disk with 
a hot spot (Abramowicz 1990). 

In chapter 6 of this thesis we have shown that there is evidence from the Seyfert 
galaxy MCG -6-30-15 that the ‘hard tail’ is variable on time scales similar to that of 
the minimum light crossing. One must determine if the disk can remain stable when 
such variations are present. 

Some evidence for inhomogeneities in the central engine of AGN already exists 
from X-ray observations (see chapter 6). The variability and spectral variability 
of the X-ray emission, which accounts for most of the emission from radio quiet 
AGN, require that the emitting region be unstable. Homogeneous and isotropic (with 
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respect to the SMBH) accretion, or stable thin disk accretion, cannot be unstable 
unless perturbations are introduced. The model presented in chapters 2, 3 and 4 of 
this thesis is natural in that the source of the power law emission which requires some 
form of magnetic fields is consistent with the fields that confine inhomogeneities. The 
same fields are also responsible for shocks that can explain spectral variability, the 
creation of clouds and outward transfer of angular momentum! This outward transfer 
of angular momentum is not inconsistent with the instability of the disk, because its 
mechanism is not stationary. The mechanism for accelerating particles (Blandford 
and Eichler 1987) is naturally explained in the context of clouds and shocks. 

It is therefore evident that the transfer of angular momentum property can be 
explained with clouds rather than disk. For the case of low accretion rates, however, 
we speculate that the disk-corona model is applicable. The critical accretion rate in 
AGN with shocks is calculated in Sivron and Tsuruta 1995. 

7.1.2 Other Results for the Cloud Model 

In deriving the physical properties of the cold clouds component we have shown 
that the optical depth in the central engine is nearly one when the luminosities are 
nearly L E dd • It is therefore reasonable to assume that the clouds are indeed optically 
thick, as long as the number of clouds in the line of sight is relatively small. This was 
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shown to be a reasonable assumption when we considered shocks as a mechanism for 
*\ the creation of these clouds (see section 4.2.2). 

The radial dependence of cloud parameters found from equations (2.11), (2.13), 
(2.14) and (2.17) is T cld oc R~ l/2 , n t h oc R ~ 2 , n cW oc R _1 , and l M oc R}! 2 . Since 
these results change with the dynamics the accuracy is ±1/2 for the exponents. Note 
that the density needed for LTE closest to the central object is approximately n t h = 
? 10 16 cm -3 for a Bolometric luminosity of 10 u erg/sec expected in a typical Seyfert 

I galaxy, whereas the broad line region clouds at ~ l0 4 R g have a density of n t h ~ 
10 lo cm -3 . From this observational result we also get oc R~ 2 . When we include 
the theoretical results of Celotti, Fabian and Rees 1992 in which the number density 
of a cloud near the SMBH is approximately n t h = 10 16-18 cm -3 for an absorption 
column of Nh < 10 21 cm -2 (for a Seyfert 1 galaxy), we find that n t h cc R~ 2 - R~ 2 66 , 
In the case that very few clouds exist one expects the outcoming radiation to 
be partially polarized, because of the effects of the source-cloud-observer angle on 
this radiation. When X-ray observation of polarization are available the results will 
probably help determine the geometry of the central engine. (Polarization which is 
observed from NLR clouds and electron scattering regions that probably reprocess 
emission from the BLR region in some Seyfert 2 galaxies was used in the determina- 
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tion of the geometry in the outer regions of Seyferts by Antonucci and Miller 1985). 
Calculations of the polarization in such a case are relatively straight forward, and 
should include clouds with Ti dd < 1. 

With the same value of /, the major spectral effects of the cloud configuration, 
as compared with the slab geometry, are (i) the increase of the UV bump, and (ii) 
the decrease of the reflection features. The difference in the reflected component is 
due to the part of the reflected component which is absorbed by the clouds between 
the reflector and the observer. Using the R dependence of the physical variables 
shown earlier we find that reflection is most effective at intermediate radii. At a 
representative radius of R = 30 R g , with Rmax ~ 300R S and black hole mass M — 
1O 7 M 0 we get T dd ~ 3 x 10 5 K and n dd ~ 10 16 cm“ 3 . For the lower limit on the radial 
dependence of the clouds size and number density, and for / ~ 0.9, we obtain l dd = 
10 12 cm and n cld = 10 -39 cm -3 for NGC 4051. Most of the covered portion of the input 
spectrum is absorbed by the clouds and reradiated as thermal radiation producing 
a large UV bump. We calculated the total integrated bolometric luminosity. Using 
equation (3) it is found that for the most extensively studied rapidly variable Seyferts, 
NGC 4051 and NGC 6814, the total primary bolometric luminosity L bd is still less 
than or roughly equal to the Eddington limit. We assumed a large /, e.g. ~ 0.9, and 


M bfl ~ 10 6 10 7 Af 0 , reasonable values. Broadening of the UV bump is evident, 

with the bump extending from the UV to the soft X-ray excess region. 

The generalization of the code which was used in chapter 3 is relatively straight 
forward. The complete code includes an extended non-thermal source, that is depen- 
dent on the distance from the SMBH, and which will therefore be represented by a 
continuous emission radius dependent function multiplied by the fraction filled by the 
hot phase. This emission will be added to the macroscopic radiative transfer equation. 
This code is not applicable in this thesis because the dominant emission mechanism 
used is bound free emission and absorption which is applicable from ~ 30/2 J} just 
outside the region in which most of the non-thermal radiation is probably produced. 

The cloud model as presented in chapters 2 and 3 might also be applied to galac- 
tic x-ray sources which are black hole candidates, because the underlying physical 
mechanisms governing these objects are considered to be similar. In the case of the 
x-ray galactic source Cygnus X - 1 our cloud model agrees better with observations 
than disk models, because the weaker and symmetrical Fe line feature predicted from 
our models with / << 1 is consistent with the observed data, while the slab reflec- 
tion models predict this feature to be stronger than and different in shape from the 
observed Fe line (Done et al. 1991). 
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In reality the clouds are most likely distributed unisotropically, forming e.g. a 
quasi-spherical flow or a torus, due to the expected presence of angular momentum. 
These geometrical effects are currently under investigation (Sivron and Tsurua 1994). 
Also, because of the plasma and magnetic effects the shape of the clouds is most 
likely non-spherical, e.g. melon-seed or sheet like, or filamental. We predict that our 
results presented here are still qualitatively adequate, because random orientation of 
the clouds will average out most geometrical effects other than the ones considered 
in this paper. Generalizing the code for such cases is best done by assuming elliptical 
disks of random and magnetic field modified shapes. 

7.2 DYNAMICAL CONSIDERATIONS 

Dynamical equations for the two components cloud model may include continuous 
equations for the optically thin component, and Maxwell-Boltzman type equations 
for the optically thick component, as is customarily done in galactic mergers (see 
Barnes and Hernquist 1992). In our case, however, we do not know the shape of 
the clouds, and their possible preferred direction. It is also hard to argue that the 
cold component clouds are much smaller than the spacing between such clouds. We 
therefore solve a simple problem, and point out the method one should employ for a 


more dynamically accurate picture. 

7.2.1 Results and Conclusions for Shocks 

A star orbiting an SMBH can be a curiosity, or a rule — depending on future 
observations, but shocks are probably present in the central engines of AGN, and 
they probably modify the emission and should be taken into account in future work. 

Natural sizes of perturbations in the outer thin Keplerian disk must be much 
smaller than R, because the speed of sound in thin disks is much smaller than the 
Keplerian speed (Frank, King and Raine 1992). The thickness of the thin disk may 
be a natural size for the perturbations. From equation (2.18) the natural size for the 
perturbations is l dd , which is similar to the size of a main sequence star in Seyferts. 

If shocks are formed the following description is the result of the treatment pre- 
sented in chapter 4: Considering the velocity of the inhomogeneity found in section 
2, and the fact that the post shock matter moves at subsonic velocity, this result 
determines that the clouds form a ‘foamy’ wall of thickness ~ c,t, ~ lO 10 cm. A 
typical filling factor of C ~ 0.01 is reasonable for the ratio of ambient and thermal 
matter under the assumption that half the ambient matter becomes thermal in the 
post shock region (with the necessary covering factor ~ 1, see Lightman and White 
1988, Celotti, Fabian and Rees 1992, Sivron and Tsuruta 1993, Bond and Matsuoka 


1994, Goerge and Nandra 1994). Therefore, on the average, one sheet of thickness 
10 8 cm can occult the continuum source. Some internal absorption is therefore ex- 
pected, but there is not yet enough data to determine whether flares or absorption 

are dominant. 

We note that in the case of virial magnetic fields the pressure on the inhomogeneity 
is similar to that calculated in equation (4.25) and cooling is independent of speed. 
This results in a 'shooting star effect’: When an inhomogeneity plunges into the 
central engine it slows down while creating tubes, until it reaches lower super Alfvenic 
speeds and a shock is formed. 

We also note that the result of the treatment in section 4.2.5 leads to a prediction 
on correlation between the slope of the power spectrum and the overall luminosity, 
because the total number per second of cold inhomogeneities accreted is proportional 
to M oc L. This is shown to be observationally plausible by Green 1993. 

The idea of optically thick shock fronts may result in periodic spectral variability 
in galactic sources, such as the Cataclysmic Variable (CV) in the field of view of NGC 
6814, is found to be consistent with observations (Sivron, Tsuruta and Leighly 1994). 
We note that CVs may be fueled by accretion disks with spiral shocks (Morfill, Spruit 
and Levy, 1993). We also note that shocks are considered to be the reason for the 
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spiral structure in the outer disk of M87 (Chakrabarti, 1995). 

7.3 RESULTS AND CONCLUSION FROM OB- 
SERVATIONS 

The results of chapters 2, 3 and 4 may be related our observational results of 
chapter 6. Some of our results support the existence of a cold or warm component 
in or near the central engine of MCG -6-30-15. Some interesting variable features 
in the spectrum can be explained in this way. Compton reflection from cold matter 
and emission from hot matter with warm (partially ionized) absorbing matter model 
can be fit to the data if both the cold and warm absorbers are near the central 
engine. This scenario is consistent with both the observational reflection and leaky 
warm absorber and the theoretical disk corona and the two component cloud models. 
Partial covering of the emission from the hot matter is with very warm matter (in 
which all elements except Iron are mostly ionized), which is also successful, may be 
checked by Asca, XTE and AXAF missions. This model is consistent with the two 
component model, but not the disk model. From the results of chapter 6 it is evident 
that there is correlation between the low energy absorption feature and the hard x-ray 
tail, and that they vary together on time scales as short as hours. Some information 


on the ionization state of this matter may already be obtained with the above results. 

See, however, the discussion in appendix F. 

In future work the variation of the lower energy part should be checked with 
the long observations that was performed this summer (1994) by the Asca satellite. 
The X-ray Time Explorer (XTE) should be invaluable in determining the spectrum 
in the now unchartered 20 - lOOkeV range, which will help establish a clear differ- 
ence between reflection and absorption models, and will determine which pair model 
is reasonable for the central engine. Possible observations of polarization with the 
calorimeter on the Japanese Astro E mission may reveal evidence for a preferable 
geometrical configuration of the central engines of AGN. (Observations of highly red- 
shifted QSOs by HST may reveal some information on the spectrum in EUV, albeit 
the observed luminosity may come from anomalously strong AGN with different spec- 


trum.) 
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Appendix A 

Reasons for Super Massive Black 
Holes 

AGN are not only among the most luminous objects in nature, but also among 
the most efficient accelerators of matter and among the most dynamically stable as- 
trophysical objects. For example, some AGN emit jets of matter that are so energetic 
and highly collimated that they extend to a distance eleven orders of magnitude 
larger than the central engine and keep their directionality throughout. Radio emit- 
ting lobes, that are almost certainly the result of the interaction of the jets with inter 
galactic material, are found either at one side, or the two opposite sides of the host 
galaxy. Theses lobes may be as far as 10 billion light years away from the nucleus 
of the host galaxy. The energy contained in those lobes exceeds the total output 
which can be produced by the central engine over a period of 10 billion years. The 
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existence of these features therefore implies that the jet emitting mechanism keeps 

its directionality for at least 10 billion years. 

The release of gravitational energy may also be possible by the extraction of energy 
directly from the ergosphere of a rotating Kerr SMBH. In the ergosphere particles 
cannot remain at one point because of the gravitational pull of the rotating black 
hole. Particles which manage to escape this region gain energy, which reduces the 
rotational energy of the black hole This process is called the Penrose process. In a 
variation on this process an elect romagneti cal torque is applied on the rotating black 
holes ergosphere (Blandford 1990). The result of such processes may be the outstream 
of energy in jets and winds. 


Appendix B 
Radio Loud AGN 


Thermal particles in the ambient gas can be accelerated to high velocities by 
shock-Fermi processes (Blandford and Eichler 1987) in the shock front. The distri- 
bution of Lorentz factors for these particles is of the form /(7 l) oe 1l + ^ * ^ or 

non-relativistic shocks, where r ~ VjV d , the shock compression, is of order 4, and 
V u and V d are the upstream and downstream (post-shock) velocities. (Relativistic 
shocks steepen the exponent even more.) Particles in the central engine are therefore 
efficiently accelerated across the shock (reference from 1995 ). 

The bremsstrahlung timescale for cooling with no magnetic fields at all is therefore 
too large to account for post-shock matter cooling. But note that some magnetic fields 
must exist in the central engine, because the accretion disk is not dense enough to 
allow for the diffusion of the ambient stellar magnetic fields away from the accreting 
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matter. 

Magnetic fields change the problem by modifying both the post-shock radiation 
and the thickness of shocks. Considering strong hydro-magnetic shocks The radia- 
tive cooling is always shorter for the compactness parameters obtained from table 
1. The compactness in that region is increased by an order of magnitude, and the 
timescales for pair-creation are very short. Matter is probably cooled in the post- 
Q shock region down to LTE temperatures. In the rarefaction region behind the shock 

this ‘cold’ matter is probably sheared into clouds by strong turbulence, and may even 
be magnetically confined. We note that collisionless shocks are much narrower, and 
therefore the time which the post-shock matter spends in the compressed region is 

O fry 

much shorter. The cooling timescale for synchrotron radiation, t coo i °c B , may 
r be decreased by an order of magnitude because of an order of magnitude increase in 

the post-shock magnetic fields, (such increase is possible in strong shocks. See, for 
example, Eilek and Hughes 1992). The synchrotron radiation alone is probably not 

c 

enough to account for the radiation. Compton cooling should be much faster than 
synchrotron. The mean free path for collisions is reduced because of the increased 
^ post shock densities and synchrotron radiation density U a . Typical Compton-cooling 

timescales, t c oc l/(nU,), are therefore reduced by up to two orders of magnitude. 
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The free-free cooling in the post-shock matter is also increased by up to two orders of 
magnitudes (see equation ( 2 ....))- The size of the post-shock matter should therefore 
be at least 10 ~ 2 R X in order that the shock emission be similar to the total source 
emission. The timescale t h for heating up the post-shock matter was derived in section 
4.1.3. The size of the post-shock matter region 

Rp» ~ v+th (BT) 

is therefore bounded by t a >t h > t c . It is apparent that there is a still a discrepancy 
here. The shock itself is the source of the power-law high energy particles. This only 
reduces the timescale by one exponential power. We therefore invoke a scenario in 
which the increased radiation results in the collapse of the hot post-shock matter into 
cool LTE matter, as in chapter 2. 


Appendix C 

The Eddington Approximation for 
an Opaque Slab 


We use the two stream Eddington approximation, as is hinted in LW88, and follow 
the method outlined in Rybiky and Lightman 1979. 

Assume the picture presented in figure 2 the radiative transfer equation in the 

case of an optically thick slab is solved in the two regions: s < 0 in which the source 
function and absorption are zero, and s > 0 in which the source function may be 
written as S v = Changing variables ds = dr/(/i(a v + <7„)), where [i = cos6 

is cosine the angle between the incoming radiation and the plane at s = 0, and s v 
is the cross section, and integrating over -1 < /x < 1 the radiative transfer equation 

becomes: 

1 dE / . = J V - S„, (C.l) 

CL V + <7„ dz 
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where J v = (1/2 ) J I v dp, H „ = (1/2) and the radiation pressure is K v = 

(1/3) J v = (1/2) f n 2 I v dn. Manipulating those equations one can easily get a second 


order equation for J u : 


<PJv =J U -S V 


(C.2) 


3(a„ + cr u ) dz 2 

In a random walk process an average single scattering albedo is 1 — e„ = <y v f{o. v + 


s„), and we get: 


44 = 3 e(J-B). (C.3) 

dr 1 

The general solution of the above equation is J = C\t^ tT + C' 2 e _v ^ T + B. Boundary 
conditions must be provided now. 


The Eddington two-stream approximation follows: the equation can now be solved 
by assuming that the entire radiation fields can be represented by radiation traveling 
in two angles, and may be written as /„ = l/y/$) + + 1/>/3)> The mean 

intensity, flux and radiation pressure are then J = / + +/ - , H = l/(2\/3 )(I + + I ) and 
I< = (l/6)(/ + + I~) = (1/3)J. The Eddington two-stream assumption is justifyable 

if there is small dependence on the angle for small angles. 

Using the relations between J, H and K we find 


J+ = J + 


1 dJ 


y/^dtau' 


(C.4) 


214 


and 


i~ = j- 


1 dJ 


(C.5) 


y/3dtau' 

The boundary conditions can now be used: I~(t = 0) = 5, n and, with the fact that 
the mean intensity cannot explode at oo but t >> 1 for an optically thick slab, the 
first term in the general solution vanishes. We thus find Ca = e~'^ T (Si n -B)/(l + y/e), 


and 


J-B = — — 


1 + V~e 

and plugging into equation (C.4) (for / + ) we get 


(C.6) 


7+ = 2B 


2 + y/e 

1 + \/c 


+ S in 


1-yS 
1 + ^' 


(C.7) 


The first term which includes the thermal emission can be treated separately, as in 
equation (3.11) in chapter 3, when the full geometry of the cloud is given, and we are 
thus left with the result given in thfe paragraph following equation (3.10). 


Appendix D 

The Case of Dense Clouds Inside 
a Disk 


In figure 3-8 we see a configuration which includes both the disk structure, which 
may come about due to the the disruption of a thin disk due to thermal instabilities 
(see chapters §1 and §4). This exercise is done in order to see how the opening angle 
affects the outcoming spectrum. 

The angle between the observer, the reflecting cloud and the point source is found 
from the product of the cloud location vector (r sin 6 cos <£, r sin 6 sin <f>, r cos 6) and 
the direction of the observer (sin /?,0,cos/?). This is one of the modifications needed 

in the code. 
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Appendix E 

The Method of Moon Phases 


The fraction of the incident luminosity which is reflected is assumed to be the frac- 
tion of the area which is observed in this case. In figure 2-..., using l = yfx 1 + y 2 
and the tan# = zjy we find that the terminator s equation is 

y — y /^/2 _ x 2 )/(i _|_ tan 2 6) = cos 6\J (l 2 — x 2 ), and the total observed area in the x , y 
plane is thus 

A = \+2[ l r° a8{l X) dydx = \nP{ 1 + cos#). (E.l) 

J i=0 Jy = 0 ^ 

Another approximation of the reflection uses the total reflected area, rather than the 
cross section: 

A = r r 6 I 2 sin 6'd$'d<j> = jt/ 2 (cos# + 1). (E.2) 

J<t >= 0 J$'=0 

The fraction of the maximum observable area is the same, but this result can be easily 
generalized to the case of ellipsoids which can roughly represent any three dimensional 
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body by varying the three semiaxes. 


Appendix F 

The Treatment of a ’Warm 
Absorber’ 


By requiring that radiative equilibrium is reached in ionized matter we get from 
the rate equation: 

jV r+ l Tl e (p) [j4 r+ i,r + u(l/ r , r +l)#r+I,r] = N T u(l/ Ti r+l)B r ,T+l (F-l) 

the Saha equation for N r+ in e (P)/N r . 

In our work we consider a simplified situation, in which the only two states involved 
are the completely ionized atom, and an atom with one electron, which is in the 
ground state as is proven in common texts (see Rybiki and Lightman). In such a case 
the induced and absorption Einstein coefficients are related through B TtT +i/ B r +i,r — 
g T+l g a ATrp 2 /g T . The induced and spontaneous coefficients are related through the 

quantum mechanical relation (?) A r +i, r = 2hi/ 3 B r +i ir /c . 
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We next consider the radiation energy density u(i/ r , r+ i) = V c where ^ = L */( R ) 
is the specific intensity of an isotropically emitting source at a distance R. Assum- 
ing the material is illuminated by a power-law source which dominates the thermal 
emission (at least as far as the ionized species are concerned) we get 

L = p ' Av-'dv (F.2) 

* Vmin 

which makes it able to explicitly calculate A = - I'min)- Since the distance 

from the central source is hard to determine it is appropriate to switch to another 
variable, the ionization parameter S = L/(4irR 2 n e chv T , r+l ). The radiation energy 
density is thus, for example, it = En e h/Q = c.D for 7 = 1, where Q = \n{u max ju mtn )- 
Plugging into the saha equation one gets 

i m\l/2 r? /*lmkT 

(F.3) 


1 V r+1 n, _ 

. Nr 9r 1+ 2l/ 3 (?/(c 2 Hn e ) 


In table .... we notice that the equivalent ‘cold’ hydrogen absorption column in the 

~ 0.8keV range is 10 - 100 times smaller than that in the 7 - 9keV region. We use 

equation (6.7) above for both of these regions, and take the ratio: 

Nr ^OWAO) 1 + jfxKgfVg (F.4) 

JV r+1 (Fe)/iV r (Fej ~ 1 + Wx 3 i(0)/Z ' 

Note that the above example is a very special case, which should be recalculated when 


competing cooling processes are present. 


Previous authors have used software such as ‘CLOUDY’ (Ferland and Rees 1988 
and references therein) for the case of MCG -6-30-15, and found absorption column 
of order ~ 10 22 cm -2 for the warm absorber and ionization parameter of order — ~ 40 
(see Nandra 1989, Fabian et.al. 1994). CLOUDY, which was extensively utilized for 
the BLR, may not be reliable for processes which occur very close to the central engine. 
In particular the above program definitely uses the density of the BLR implicitly when 
calculating the ionization parameter. Higher densities will result in different cooling 
and heating mechanisms. Overfitting of emission/absorption features which are not 
efficiently detected by Ginga may also be a problem with such programs, as they 
tend to distort the spectral features which are actually observed. Software under 
development for the region (Netzer 1993) may be utilized for our subsequent work. 


